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Abstract 

cn 

^H We construct a two-dimensional crystal melting model which reproduces the 

>► BPS index of D2-D0 states bound to a non-compact D4-brane on an arbitrary 

•i-H 

^ toric Calabi-Yau singularity. The crystalline structure depends on the toric 

divisor wrapped by the D4-brane. The molten crystals are in one-to-one corre- 
spondence with the torus fixed points of the moduli space of the quiver gauge 
theory on D-branes. The F- and D-term constraints of the gauge theory are 
regarded as a generalization of the ADHM constraints on instantons. We also 
show in several examples that our model is consistent with the wall-crossing 
formula for the BPS index. 
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1 Introduction and summary 

The geometry near D-branes is probed by lighter branes bound to them. One of the 
most well-known examples is that Dp-branes bound to D(p + 4)-branes on the orbifold 
C 2 /Zat describe instantons on the resolved A^-i ALE space [I]. 



Recently, there has been remarkable progress in the study of D0-D2 states bound to 
a D6-brane on a toric Calabi-Yau three-fold. From the above viewpoint, such D-branes 
probe the Calabi-Yau geometry wrapped by the D6-brane. In fact, the BPS index of 
the D-brane bound states is evaluated by counting the molten configurations of a three- 
dimensional crystal [21 El HI El E] , whose crystalline structure is determined by the toric 
diagram of the background Calabi-Yau three-fold [5jjj Moreover, the thermodynamic 
limit of the molten crystal describes the smooth geometry of the mirror Calabi-Yau 
three-fold [HI [12]. This suggests that the melting crystal gives a "discretization" of 
the background Calabi-Yau geometry. 

The generalization of the crystal melting model to D4-D2-D0 states has partially 
been studied. In [131 HI], the authors considered D2-D0 states bound to a non-compact 
D4-brane on a divisor of the (generalized) conifold, and constructed a two-dimensional 
statistical model which reproduces the BPS index of the D4-D2-D0 states. After the 
success of the D6-D2-D0 crystal melting, it is natural to expect that the structure of 
the two-dimensional model is related to some property of the toric divisor wrapped by 
the D4-brane. However, such a relation has not yet been clarified. The main reason 
for this is that the prescription given in [T31 E] is ad hoc and not derived from the 
BPS condition for the D-brane bound states. 

In this paper, we derive a general method to construct a two-dimensional crystal 
melting model for D4-D2-D0 states on an arbitrary toric Calabi-Yau three-fold, by 
solving the BPS condition for the D-branes. We put a D4-brane on a non-compact 
toric divisor P of a toric Calabi-Yau three-fold Y, and count BPS D2-D0 states bound 
to it. Here the D2-branes are wrapped on compact two-cycles of Y, and the DO-branes 
are point-like in Y. We particularly consider the singular limit of Y, in which D2-D0 
states are realized as fractional branes localizing at the singularity. We identify the su- 
persymmetric gauge theory on the D-branes, and solve the F- and D-term constraints. 
There is a natural torus action on the moduli space A4ba of supersymmetric vacua, 
which essentially comes from the toric action on Y. We then show that the torus fixed 
points of .Md4 axe in one-to-one correspondence with the molten configurations of a 
two-dimensional crystal. The crystalline structure depends on the choice of the divisor 
wrapped by the D4-brane. 

To give a short summary of this paper, let us first consider the simplest case of 
Y = C 3 and T> = C 2 C Y. Since C 3 has no compact two-cycle, we can only consider 

There are also several works on the crystal melting description of the wall-crossing phenomena 
[TUBUS] and refinement [ID]. 



Figure 1: A Young diagram (left) can be regarded as a molten configuration of a two- 
dimensional crystal (right). The crystal infinitely extends in the upper-right region. 



D4-D0 states. The low-energy effective theory on k DO-branes bound to a D4-brane 
on T> is a d = 1 supersymmetric U(k) gauge theory with 8 supercharges. In the 
d = 4,J\f = 2 language, the theory includes three adjoint chiral multiplets B a for 
a = 1, 2, 3 which come from DO-DO strings with the usual M = 4 superpotential 



ti(B 1 [B 2 ,B 3 \). 
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There are also a fundamental and an anti-fundamental chiral multiplet /, J which 
come from D4-D0 strings with the J\f = 2 superpotentiajj 
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The well-known fact is that the F-term conditions for the above superpotential imply 
the ADHM constraints on /c-instantonsJ3 

[B 1 ,B 2 ] + IJ = 0. (1.3) 

Here, the effect of the non-compact D4-brane clearly appears in the additional super- 



potential (1.2). The moduli space A^D4 of supersymmetric vacua admits a natural 



torus action. The BPS index is essentially equivalent to the number of the torus fixed 
points of A4d4, which are labeled by Young diagrams [151 HB]- Note that each Young 
diagram is regarded as a molten configuration of a two-dimensional crystal composed 
of square boxes (figure [l]). Therefore, counting the torus fixed points is equivalent to 
counting the molten crystals. The aim of this paper is to generalize this to arbitrary 
Y and V. 

When Y is a general toric Calabi-Yau three-fold and T> is one of its toric divisors, 
the low-energy effective theory on D-branes is identified by using the technique of 



5 Here we assume without loss of generality that B3 describes the fluctuations of the DO-branes in 
directions transverse to the D4-brane. 
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The real part of the ADHM constraints comes from the D-term constraint. 



brane tiling [S]. In section [21 we briefly review the work [5j on the crystal melting 
model for D6-D2-D0 states on a toric Calabi-Yau three-fold, including the technique 
of brane tiling. The brane tiling was originally developed in the study of the D-brane 
construction of J\f = 1 supersymmetric gauge theories (TTJ [181 [191 [20] (F° r reviews, see 
also J2TJ 122])- We particularly study brane tilings with so-called isoradial embedding, 
as mentioned in subsection 13.11 The main difference from the D6-D2-D0 case is that 



the D4-brane induces an additional superpotential such as (1.2). In subsection 3.2, we 
use a technique developed in (23] to identify such an additional potential for arbitrary 
Y and T>. We then claim in subsection |3.3 that the F- and D-term constraints imply 



that some chiral multiplets have vanishing vev's on supersymmetric vacua. In the 
above example of Y = C 3 , B3 is such a chiral multiplet. The set of chiral multiplets 
with vanishing vev's is specified by a so-called "perfect matching." The proof of our 



claim is given in subsections 3.7 and 3.8, where we use the results of [24j. Based 



on this observation, we show that the resulting moduli space A^D4 is embedded in 
the moduli space A^D6 of the parent D6-D2-D0 state. Here, the parent D6-D2-D0 
state is obtained by replacing the D4-brane with a D6-brane wrapping on the whole 
Calabi-Yau Y. The inclusion map i : A^D4 *->■ -Md6 is characterized by the perfect 



matching. We then show in subsection 3.4 that the torus fixed points of A^D4 are in 



one-to-one correspondence with the molten configurations of a two-dimensional crystal. 
This two-dimensional crystal is, in fact, a "slope face" of the three-dimensional crystal 
associated with the parent D6-D2-D0 counting. Moreover, choosing a different divisor 
T> of the same Calabi-Yau Y gives a different slope face of the same three-dimensional 



crystal. As described in |3.5[ the boundary of the two-dimensional crystal is given by 
so-called "zig-zag paths." 

In section [4], we give several examples in which Y is C 3 , the conifold, the suspended 
pinch point, and the orbifold C 2 /Zjy x C In particular, when Y is the conifold, our 
crystal melting model reproduces the triangular partition model proposed in [13]. If 
Y is the suspended pinch point, our model reproduces the oblique partition model 
proposed in [T4] . Furthermore, when Y = C 2 /1*n x C and T> = C 2 /Zjv C Y, our 
model reproduces the orbifold partition model [251 ES, [27] whose partition function 
agrees with the M = 4 U(l) instanton partition function on the An-i ALE space, 
i.e. the level-one character of the affine SU(N) algebra. Note that this example is 
the original setup of [T], which has been mentioned at the very beginning. In fact, 
our setup is a generalization of that of [1] to an arbitrary toric divisor. In some of 
the examples, we explicitly show that our model is consistent with the wall-crossing 



formula for the BPS index. 

One of interesting future works is to extend our result to the multiple D4-branes 
with D2- and DO-branes. When Y = C 2 /Z N x C and V = C 2 /Z N C Y, this setup 
produces the affine Lie algebra character with a higher level than one. It will also be an 
interesting future problem to assign certain weights to each fixed point of U(l) 2 in the 
moduli space and study an analogue of the Nekrasov's partition function (T6J [28] . This 
generalization should have a good realization in the string theory and the M-theory. 
Furthermore it should be related to some observables of a 2-dimensional conformal 
field theory by the AGT relation [2H] . From the viewpoint of the gauge theory on the 
D4-brane, this work can be thought of as the instanton counting in the d = 4, M = 4 
supersymmetric gauge theory on a toric divisor|_] There are several interesting works 
[3TI l32l 133] on the crystal melting in this context. It is worth studying the relation to 
these works. It would also be interesting to apply our method to a Calabi-Yau three- 
fold with compact four-cycles. Although we do not explicitly consider such an example 
in this paper, the application is straightforward. When Y contains compact four-cycles, 
the melting crystals also count the charge for compact D4-branes. In particular, it 
would be interesting to study the relation to the work of [31]. Another interesting 
direction would be to study the connection to the works [351 EE1 E3 EHJ [39] , where the 
crystal melting model was studied in the context of five-dimensional supersymmetric 
gauge theories. It is also worth studying the relation to the works [4"0l I4"lj |4*2| |4"3] on the 
BPS chiral operators of gauge theories on D-branes at toric Calabi-Yau singularities. 

2 Crystal melting for D6-D2-D0 states 

We here briefly review the crystal melting model for D6-D2-D0 states on a toric Calabi- 
Yau singularity, mainly following [5]. We consider the BPS bound states of a single 
non-compact D6-brane and arbitrary numbers of D2 and DO branes, where the D6- 
brane is wrapping the whole Calabi-Yau, the D2-branes are wrapping compact two- 
cycles and the DO-branes are point-like in the Calabi-Yau three-fold. 



7 To be more precise, the theory is a topologically twisted Af = 4 super Yang-Mills theory. The 
fact that the divisor is embedded in a Calabi-Yau three-fold implies that the topological twist is the 
Vafa-Witten type twist [50] . 



"1 

<> — I — n 

4 2 

* ► 

O 1 O 

-3 




Figure 2: Left: The toric diagram of the conifold. Each line segment s corresponds to 
a non-compact curve j3 s in the conifold. Right: The corresponding T 2 -fiber. We have 
a non-degenerate T 2 -fiber at a generic point of the toric base 1R 3 . However, the T 2 
degenerates to S 1 along some semi- infinite lines from p^. The semi-infinite lines are 
in one-to-one correspondence with the curves j3 s . At the intersection p%, the T 2 -fiber 
shrinks into a point, giving rise to a singularity. 



2.1 Quivers on D2-D0 from brane tilings 

A toric Calabi-Yau three- fold Yy, is roughly regarded as a (T 2 x IR)-bundle over M. 3 , 
where the T 2 -fiber degenerates in a subspace specified by a toric diagram S. We 
mainly consider Yg without compact 4-cycles. The toric diagram E is a convex lattice 
polygon, in which every vertex is associated with a toric divisor of Yj> Every line 
segment s in £ is associated with a non-compact curve (3 8 in Y"s, along which the 
T 2 -fiber degenerates to S 1 . The degenerate cycle is specified by the slope of the line 
segment s. Namely, if s is stretched between two vertices (pi,qi) and (p 2 , 92) in S, 
then (qi — q 2 , —p\ + p2)-cycle of T 2 degenerates along the curve f3 s . The transverse 
(j»i — P2,qi — q2)-cjcle then generates an isometry of the curve (3 S . At the intersection 
ps of all the curves j3 s , the T 2 -fiber shrinks into a point, giving rise to a singularity 
(See figure [2]). 

Let us consider DO-branes localized at the singularity p-%- In order to evaluate the 
supersymmetric index for such BPS DO-branes, we identify the low energy effective 
theory on the DO-branes. We first take the T-duality transformation along the two 
directions of the T 2 -fiber, which maps the DO-branes to D2-branes wrapping on the 
whole T 2 . On the other hand, the toric Calabi-Yau geometry itself is mapped to inter- 
secting NS5-branes in flat spacetime |HJ HHJ HHJ S6] because the T-duality exchanges 





Figure 3: Left: Each NS5-brane is a semi-infinite line in the toric base and wrapped 
on a one-cycle of the T 2 -fiber. Right: All the NS5-branes intersect at the point p^, 
where we have a brane tiling system. The red lines are the NS5-branes, which "bend" 
to fill up the gray regions. In the white regions, we only have dynamical D-branes 
without NS5-branes. Each blue arrow is associated with a massless chiral multiplet. 



the source of the KK gauge field with that of the NSNS B-field. To be more concrete, 
the curve /3 S for every line segment s of E is mapped to a single NS5-brane wrapped 
on a semi-infinite tube. The tube is an S^-fibration over a semi-infinite line from ps 
in the toric base. If s is stretched between (pi,q±) and (^2,^2), then the S 1 is the 
(pi — P2i Qi — <?2)-cycle of the T 2 -fiber (See figure [3]). The NS5-branes also extend in 
M 4 transverse to the six-dimensional space we are considering. 

The D2-branes wrapping on T 2 , which come from the original DO-branes, are now 
divided into several "tiles" by the intersecting NS5-branes (figure [3]). Here, the conser- 
vation of the NS5-charge implies jlTj that the intersecting NS5-branes in fact "bend" 
to fill up some of the tiles (The gray regions in figure [3]). To be more precise, if the 
boundary NS5-branes of a tile gives a definite orientation, the tile is filled with an 
NS5-brane. Such tiles are classified into two types, depending on the orientations of 
the NS5-brane. In figure [3j the dark and light gray regions are filled with NS5-branes 
with opposite orientations. The white regions has no NS5-branes, and are filled only 
with the D-branes. By construction, neither two gray regions nor two white regions 
share any edge. Furthermore, the corner of a light gray region is attached to the corner 
of a dark gray region, and vice versa. 

Now, we consider the low-energy effective theory on the D-branes. When we reduce 
the two directions of T 2 , we obtain a supersymmetric quantum mechanics on the 
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Figure 4: Left: The quiver diagram Q associated with the brane tiling in figure [3| The 
diagram is defined on T 2 ; The upper (right) and lower (left) dotted lines are identified. 
Here Nq and N 2 are the DO- and D2-charges on the conifold, respectively. Right: The 
left quiver diagram is equivalent to the well-known quiver for the conifold. 



world-volume of the D-branes, whose field content can be read off from the brane 
configuration in T 2 . Each white region in T 2 gives a U(N ) gauge group, where N is 
the original DO-brane charge. The gray regions give rise to no gauge multiplet because 
they are filled with NS5-branes. We also have a bifundamental (or adjoint) multiplet 
at each intersection point of the white regions. Such a bifundamental is expressed as 
an arrow from one white region to the other (figure [3]). The orientation of the arrow is 
determined by the relative positions of two adjacent gray regions. We determine it so 
that the arrow goes from bottom to top when the adjacent light gray region is on its 
left sidejj This definite orientation means that the bifundamental (or adjoint) is not a 
hyper multiplet but a chiral multiplet. Thus, the low-energy theory on the D-branes 
is a d = l,J\f — 4 quiver quantum mechanics. For example, the quiver diagram for 
the brane configuration in figure [3] is shown in figure |4| For a given quiver diagram Q 
on T 2 , we denote by Qo the set of nodes in Q, as well as by Q\ the set of arrows in 
Q. The set of faces in Q is denoted by Q 2 - We sometimes denote by Q the universal 
cover of the quiver diagram Q. 

The superpotential for the chiral multiplets comes from string disk amplitudes 
associated with the gray regions. If a light gray region is surrounded by the chain of 
arrows Xi, ■ ■ ■ , X n , then there is a superpotential term of the form 



ti{X 1 ---X n ). 



(2.1) 



? Of course, the adjacent dark gray region is on the right side of the arrow. 



Here we used the same symbol X a to denote a chiral multiplet associated with the 
arrow X a . On the other hand, a dark gray region contributes 

-tr(Xi-..X m ) (2.2) 

to the superpotential if it is surrounded by the chain of the arrows X\, ■ ■ ■ , X m . The 
minus sign here is due to the opposite orientation of the NS5-branes. Since each arrow 
is attached to one light gray and one dark gray region, every chiral multiplet appears 
twice in the superpotentialjj 

So far we have only considered DO-branes at the original Calabi-Yau three- fold Ys, 
which leads to the same rank of the gauge groups. Now, let us consider additional D2- 
branes wrapping on some compact two-cycles of Y%. In the singular Calabi-Yau limit, 
all the compact cycles are vanishing and the D2-charge is realized as the fractional 
DO-charge. For example, when we originally have N 2 D2-branes on a vanishing two- 
cycle, the T-duality maps them to D2-branes which fills one of the white regions in 
T 2 , increasing the rank of the corresponding gauge group by N 2 (figure Gj). Thus, 
introducing D2-charges changes the rank of the gauge groups. 

We have seen that the T-duality maps the toric Calabi-Yau geometry to the brane 
tiling system, from which we can read off the quiver diagram Q on T 2 and the superpo- 
tential. Let us here mention the dual diagram Q y of Q. Namely, we consider a graph 
Q y on T 2 such that Q y ~ Q 2 , Q y — Qi and Q 2 ~ Q . In Q y , every face is associated 
with a gauge group while every vertex gives a superpotential term. Every line segment 
is associated with a chiral multiplet. We denote this dual map by ip : Q — > Q y . We 
sometimes denote the universal cover of Q y by Q y . We assign the color of white or 
black to a vertex / G Q v , depending on the orientation of the corresponding face 
tp~ 1 (f). If ip~ 1 (f) is a light gray (dark gray) face, then we assign white (black) to /. 
Then, it follows from our construction that the dual diagram Q y is a bipartite graph 
where any line segment connects one black and one white vertex (figure [5]). In mathe- 
matics, such Q y defines a dimer model on T 2 . Therefore, we sometimes call vertices in 
Qq "dimer vertices" while those in Qq "quiver vertices." This dimer model Q y plays 
an important role in the main part of this paper. 



3 This implies that we can set all the coupling constant in W to be +1 or —1 by rescaling the 



superficial X; L . In (2.1 1 and (2.2 1, we have already taken into account such rescalings. 
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Figure 5: The quiver diagram Q (left) and the dual graph Q w (right) on T 2 . Each 
vertex in Q y is black or white colored, depending on the orientation of the dual face 
in Q. By construction, Q y is always a bipartite graph, giving a dimer model on T 2 . 

2.2 Adding a flavor D6-node 

We now put an additional D6-brane wrapping on the whole Calabi-Yau Y^. Since 
such a D6-brane is non-compact, we can regard it as a flavor brane. In the quiver 
language, it adds an additional flavor node * to the quiver diagram Q. Note that the 
T-duality maps the D6-brane to a non-compact D4-brane localized at a point in T 2 . 
We assume that the point is attached to a white tile in the brane tiling. Then the quiver 
has an arrow I from * to an another node i G Qo, which describes a "quark" in the 
fundamental representation of the gauge group associated with i. Thus, adding a single 
D6-brane implies an extended quiver Q with Q = Q U {*} and Q 1 — Q 1 U {/}. Since 
there is no gauge invariant operators involving J, we have no additional superpotential 
induced by the D6-brane. From the mathematical viewpoint, the flavor node * gives 
a "framing" for quiver representations. 

2.3 Moduli space and ^-stable modules 

The D-brane bound states we are considering can be seen as BPS particles in M 4 , 
and our aim is to evaluate the BPS index of such BPS particles. Since the D-brane 
world- volume itself breaks half the supersymmetry, the BPS index is regarded as the 
Witten index of the quiver quantum mechanics on the D-branes. To evaluate the 
Witten index, we first have to identify the moduli space of supersymmetric vacua. 

The moduli space of the world-volume theory is parameterized by the supersym- 
metric configurations of scalar fields. Some of the scalars correspond to fluctuations of 
the D-branes in the Calabi-Yau three- fold Y"s, while the others express fluctuations in 
IR 3 . Here the latter describes the position of the BPS particle in M. 3 and we fix it. Then 
the remaining moduli space is exactly the same as the moduli space of a d = 4, M = 1 



11 



quiver gauge theory with the same quiver diagram Q and the superpotential W . Below, 
we describe this moduli space in terms of quiver representations. 

We first define the so-called "path-algebra" CQ which is generated by paths in the 
quiver diagram Q. Any path in Q starts at a node, follows some arrows and terminates 
at a node. The product of two paths is defined by connecting the tail of the first path 
with the head of the second path. Here, if the tail of the first and the head of the 
second are not attached to the same quiver node, the product is defined to be zero. 
For each node £ g Qq, there is a special element e^ of CQ which corresponds to a path 
from £ to £ with zero length. Since e^ is regarded as the projection onto the set of 
paths terminating at £, the element Yleedo ee is the multiplicative identity of the path 
algebra CQ. Physically, a CQ-module corresponds to a configuration of chiral fields 
in the quiver quantum mechanics, which might break supersymmetry. 

In order to impose the F-flat condition, we consider a quotient 

A = CQ/F (2.3) 

where J 7 is the ideal of CQ generated by all the derivatives dW/dX a for X a g Qi| 10 | 
Physically, an A-module expresses a F-flat configuration of the chiral fields. Note here 
that an A-module M has a natural grading M = (B ee n Mi where Mi = e^M. Since 
A includes a subalgebra isomorphic to C, each module M*. is naturally a complex 
vector space. For each k g Q , dimM fc is identified with the rank of the gauge group 
associated with the node k. On the other hand, dim M* expresses the rank of the flavor 
symmetry group associated with the D6-node. Since we only have a single D6-brane, 
we set dim M* = 1 . 

On the other hand, the D-flatness condition is known to be equivalent to the 9- 
stability [48J defined as follows. For a given set of real parameters 9k for all k g Qq 
and #*, the slope function of an A-module M is defined by 

9(M) = J2 OkdimM k + 0*dimM*. (2.4) 

fceQo 

Now, for a given A-module M, we fix 9^,9* so that 9(M) = 0. Then, M is called 
9 -stable if every non-zero proper sub-module M' C M satisfies 9(M') < Op] In [4~8] . 
it was shown that ^-stable modules with complexified gauge groups are in one-to-one 

10 Here we denote a path corresponding to an arrow X t by the same symbol Xj. 
11 We here use the opposite inequality sign in comparison to |48j . which is just a matter of conven- 
tion. 
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correspondence with the D-flat configurations, where 9 k, 9* are identified with the FI 
parameters of the quiver gauge theory. The condition 9(M) = is easily understood 
when identifying 9k, 9* with the FI parameters. In fact, the D-term conditions are 
schematically written as 

where Si and Tg are the sets of arrows in Q\ which start and end at I G Qo respectively. 



Summing up the trace of (2.5) for all £ G Qo, we obtain 9(M) = 0. This implies that 
there are \Qq\ independent ^-parameters for a given A-module. In this paper, we 
always set 9 so that 9k < for all k G Qo- Note that this implies #* > 0. Changing the 
^-parameters generically gives rise to the wall-crossing phenomena of the BPS states. 
For more on the stability condition, see appendix \A\ 

From this argument, we can see that ^-stable A-modules are in one-to-one cor- 
respondence with the supersymmetric vacua of the quiver quantum mechanics. In 
the next subsection, we use this correspondence to study the torus fixed points of the 
moduli space. In particular, we use the fact that any 9-stable A-module M with 9k < 
for all k £ Qo is a cyclic module generated by an element m G M*, that is, M = Am 
[7]. In fact, if Am ^ M then Am is a proper submodule of M with 9(Am) > 0, which 
contradicts with the ^-stability of M. Note that, since M* is one-dimensional, such an 
element m is essentially unique. 

2.4 BPS index and molten crystals 

We now come to the main part of this section. We evaluate the Witten index in terms 
of the ^-stable A-modules. First of all, the Witten index can be calculated via the 
localization with respect to U(l) 3 = U(l) 2 x U(1)r action on the moduli space. Here, 
U(l) 2 comes from the toric actions on the Calabi-Yau three-fold Is while U(1)r is the 
R-symmetry of the theory^] 

To see the £/(l) 3 -actions explicitly, let us consider a map t:Q\—¥ U{1). Such a map 
assigns a global [/(l)-action to each chiral multiplet. We require that t keeps all the F- 
term conditions, or equivalently, t keeps the superpotential up to an overall rescaling. 
We denote the set of all such t by T. Since the D-term constraints are obviously 
invariant under t G T, the moduli space of supersymmetric vacua is symmetric under 



12 To be more precise, £/(1),r is the R-symmetry of the parent d = 4, Af = 1 supersymmetry of our 
d = l,Af = 4 quiver quantum mechanics. 
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the action of T. Let us consider how many independent t-actions there are. For each 
face / e Q2, we define 

n, = J] t(X), (2.6) 

xedf 

where the product is taken over chiral fields surrounding /. Preserving the F-term 

conditions is equivalent to requiring that the ratio nf r /nf 2 is precisely invariant for any 

/l; fi £ Qi- Note that the ratios are not all independent because there is an identity 

J] nf n{f) = 1, (2.7) 

/6Q2 

where sign(/) = ±1 depending on the color (light or dark gray) associated with the 
face /. This identity follows from the fact that any chiral field is attached to one 
dark and one light gray region. Due to this, there are only (IQ2I — 2) independent 
ratios nfa/rifa. Requiring all of them invariant imposes (IQ2I — 2) constraints on the 
t-actions0 We then find T ~ l7(l)l<W-(l<fcl- 2 ). 

Note that some of the t-actions are absorbed into gauge transformations. In the 
U(lyQo\ subgroup of the gauge group, the diagonal U(l) keeps all the chiral fields 
invariant. The other [/(l)^ ! -1 can absorb (\Qo\ — 1) degrees of freedom of T. The 
independent global symmetry of the theory is thus T/L^l)'^ ' -1 ~ U(l) 3 , where we 
used the Euler formula for T 2 : \Q \ — \Qi\ + \Q 2 \ = 0. This means that each t- 
action is an [7(l) 3 -action modulo gauge transformations. Since ^-stable A-modules 
are generated by a single element m, this ?7(l) 3 -action is naturally extended to the 
modules, where we set m to be invariant under T. 

The £7(l) 3 -fixed points of the moduli space are equivalent to ?7(l) 3 -invariant 9- 
stable A-modules. It was shown in [1] that the latter is in one-to-one correspondence 
with so-called "finite ideals" defined as follows. We first define A* to be the set of the 
F-term equivalence classes of paths starting at the node *. For any [x], [y] G A*, we 
write [x] < [y] if there is a path z and representatives x,y of [x], [y] so that y = zx. 
Then, a finite ideal n of A* is defined as a subset of A* with the following property: 

For any [x], [y] G A* satisfying [x] < [y], if [y] G n then [x] En. 

For a given finite ideal it of A*, an A-module spanned by all the elements in n is a U(l) 3 - 
invariant ^-stable module. On the other hand, any [7(l) 3 -invariant ^-stable module 



i:s 



As we will see later, there are equal numbers of faces with sign(/) = +1 and those with 



sign(/) = — 1. This guarantees that the identity (2.7| is consistent with the overall rescaling of 
the superpotential. 
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Figure 6: Left: The periodic quiver Q for the conifold case, which is the universal cover 
of the quiver diagram Q on T 2 . Right: The three-dimensional crystal Ca, is obtained 
by putting atoms for all the elements of A*. Each finite ideal of A* corresponds to a 
molten configuration of the crystal. 

has its basis corresponding to a finite ideal of A*. This is clearly a generalization of 
the fact that the torus fixed points of the moduli space of instantons in C 2 are labeled 
by Young diagrams. 

What is important here is that the finite ideals of A* are expressed as molten 
crystals [5] . To see this, we first consider the universal cover Q of the quiver diagram 
Q, which is also attached to a D6-node * at a reference node (figure pi). We call Q 
the periodic quiver. We now start at * and follow all the F-term equivalence classes 
of paths in A*, putting an "atom" on the ending node of each path. From the dimer 
model viewpoint, we place an atom on a face of Q y . An atom on k G Qo has its "color" 
determined by p(k) where p : Q — > Q is the natural projection. When we consider two 
different equivalence classes [ci], [c 2 ] G A* ending at the same node k G Qo, we need 
to consider the "depth" of atoms. It was shown in [I] that any path from * to a node 
k G Qo is F-term equivalent to ffcUT for some £ G N, where v k is the shortest path from 
* to k and a; is a loop around some face of QJ^j Here to and v u are not unique, but the 
integer £ is uniquely determined. Thus, the elements of A* are completely classified 
by its ending node k G Qo and the integer £. Then, we determine that for a F-term 
equivalence class [v ko/] G A* we put an atom on the node k at the depth £. If we put 
atoms for all the elements of A*, we obtain a three-dimensional crystal Ca, on the 
periodic quiver Q (figure |6J). Note here that there is a bijection / : A* — > Ca„- For 

1 To be precise, this equivalence relies on some conditions on the dimer model. In this paper, 
we only consider dimer models with "isoradial embedding" as explained in subsection |3.1| For such 
dimer models, all the conditions are satisfied [11 [Ml EPJ1 150] . 
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any two atoms a,/3& Ca„, we say that there is a "bond" from a to /3 if there is an 
arrow from f~ 1 (a) to / _1 (/3). It is now clear that a finite ideal of A* corresponds to 
a sub-crystal p of Ca„ such that 

a bond from (3 G Ca„ to a G p implies /3 G p. 

Now, recall that the Witten index of the quiver quantum mechanics on D-branes 
is evaluated as a sum over the [7(l) 3 -fixed points of the moduli space. Since the fixed 
points are in one-to-one correspondence with the molten configurations of Ca„ , we can 
write the Witten index as a sum over the molten crystals p. The D2 and DO charges 
for a given p are determined by dk '■= dimM^ for k G Q$. This dk is in fact the number 
of atoms in p which are associated with the fc-th quiver node. To be more specific, let 
us define the generating function of the Witten index fl(j) as 

£bps= Y. ^(^ + m I Pi-ndV)q n l[Qf, (2.8) 

n,m'eZ / 

where q and Qj are Boltzmann weights for DO and D2 charges. The index I runs 
over 1, • • • , N where N is the number of the compact two-cycles in the Calabi-Yau 
three-fold Y^. We also define the generating function of the molten crystals as 

2 crystal = j2(-i) dimcM ™ n **> ( 2 - 9 ) 

where Xi is a Boltzmann weight for the i-th node in Q and .Md6 denotes the moduli 
space of the BPS states with charges {di}. The sign factor depends on the complex 
dimension of the moduli space -Md6- Now, what was pointed out in [5] is that 

^BPS = ^crystal (2-10) 

holds under a suitable identification between the Boltzmann weights (q,Qi) and Xj. 
The explicit identification between the Boltzmann weights depends on the original 
Calabi-Yau three- fold Yg. 



3 Crystal melting for D4-D2-D0 states 

In this section, we replace the D6-brane with a non-compact D4-brane on a toric 
divisor and consider BPS D2-D0 states bound to the D4-brane. The main difference 
is that the flavor D4-brane gives an additional superpotential term which reduces the 
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Figure 7: Left: An isoradial embedding of Q v on T 2 , shown in its universal cover. 
Every vertex / G Qq at the boundary of a face k G Q^ is on a unit circle Ck surrounding 
the face. Right: A Dimer model with an isoradial embedding implies a rhombus 
tiling of T 2 . 



moduli space of the quiver quantum mechanics on the D-branes. We will show that the 
BPS index of the D4-D2-D0 states is evaluated by counting two-dimensional melting 
crystals. We particularly concentrate on dimer models with isoradial embedding as we 



explain in subsection 3.1 



3.1 Isoradial embedding 

As reviewed in section [2j the brane tiling gives a dimer model Q v associated with 
the toric Calabi-Yau three-fold Is- The universal cover Q v of Q w plays an essential 
role in the construction of the crystal melting model. In the rest of this paper, we 
particularly consider dimer models with isoradial embedding. The isoradial embedding 
is an embedding of Q y in T 2 so that every vertex / G Qq at the boundary of a face 
k G Q2 is on a unit circle Ck (right picture of figure [7]). Dimer models with an isoradial 
embedding have been discussed in the study of R-charges in the IR superconformal 
fixed point [20] ■ The necessary and sufficient condition for Q y to admit an isoradial 
embedding is reviewed in appendix [Bj What is important here is that the existence of 
a dimer model with an isoradial embedding was shown in [51] for an arbitrary toric 



Calabi-Yau three- fold _ In the rest of this paper, we will focus on such a class of 
dimer models. 

An important property of such dimer models is the existence of a rhombus tiling 
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Figure 8: Left: After the T-dual transformations, the D4-brane becomes a flavor D2- 
brane filling up the facet £ i n the toric base, whose boundary is embedded in two 
NS5-branes along t\ and £2- Right: In the brane tiling at ps, the intersection point 
of the boundary NS5-branes is attached to two dynamical D2-branes i and j. The 
D4-node is located at the intersection, and induces two massless "quarks" / and J. 

on T 2 . Given an isoradial embedding of Q v , we can draw a line from a quiver vertex k 
to a dimer vertex / if / is on the unit circle c^. Such lines form rhombi as in the right 
picture of figure [7} Note that each rhombus is associated with a chiral multiplet X a . 
The fact that any dimer vertex is on some unit circle implies that this procedure leads 
to a rhombus tiling on T 2 . We will exploit this property heavily in the main part of 
this section. 



3.2 Flavor D4-node 

Let us now discuss the location of the flavor D4-brane in the brane tiling. We assume 
the D4-brane is wrapping on a toric divisor T> corresponding to a corner of the toric 
diagram E. The reason for this is that we want to make T> non-degenerate in the 
singular Calabi-Yau limit. The projection of T> to the toric base is a two-dimensional 
facet # bounded by two semi-infinite lines £\ and £ 2 - The lines £±,£2 are associated 
with two external legs of the toric web-diagram. The T-duality along T 2 maps the 
D4-brane to a flavor D2-brane filling up the facet $. The flavor D2-brane is bounded 
by two intersecting NS5-branes, each of which wraps on a one-cycle of T 2 and extends 
along £1 or £ 2 (figure [8]). At the intersection point p% of £\ and £2, we have a brane tiling 
system as in figure [3j Note that the boundary of the flavor D2-brane is a point in T 2 
and embedded in the boundary NS5-branes. In particular, in the brane tiling system 



at p-s, the flavor D2-brane is located at the intersection point of the two boundary 
NS5-branes (figure pi), which gives a flavor "D4-node." We use the same symbol * to 
denote the D4-node in the quiver diagram, but it gives rise to quite different physics 
from the D6-node. 

In the brane tiling, an intersection of two NS5-branes is always attached to two 
dynamical D2-branes (coming from the original D2-D0 states on Ys). This implies that 
the D4-node has two massless "quarks" attached to it. To be more specific, suppose 
that the D4-node * is adjacent to dynamical D2-branes i and j, and there is a chiral 
field Xp associated with an arrow from i to j (the right picture of figure [8]). Then we 
have a "quark" I associated with an arrow from * to i as well as an "anti-quark" J 
associated with an arrow from j to *. The quark and anti-quark are involved in the 
following additional superpotential term [23J: 

WW = JX F I, (3.1) 

which gives additional F-term conditions 

JX F = 0, X F I = 0. (3.2) 

The total superpotential is now written as W = W + Wfl avor . The first term W does 
not contain /, J, and can be read off from the brane tiling as explained in subsection 



2.1 The additional potential (3.1) also modifies the F-term condition with respect to 
Xp as 

dW 



dXi 



+ IJ = 0. (3.3) 



Thus, we now have two different F-term constraints (3.2) and (3.3) from the D6- 
D2-D0 case. Both of them are induced by the additional superpotential (3.1). In 
the next subsection, we discuss how the two differences change the moduli space of 
supersymmetric vacua. 

The location of the node i will be important in the construction of the two- 
dimensional crystal melting model. In fact, it gives a framing for the quiver rep- 
resentation. We stress here that i is the starting node of Xp. 

We also mention that the boundary NS5-branes of the facet $ in general have 
several intersections in T 2 . For example, when the Calabi-Yau is C 2 /Zjv x C and 
the D4-brane is on C 2 /Ztv, there are N different intersection points of the boundary 
NS5-branes. In general, the D4-node is located at one of the intersection points, and 
the choice is related to the holonomy of the gauge field at infinity on the D4-brane. 
We will discuss this in more detail in subsection 14.41 

19 



ra 



m 1 



m 2 



m 3 







(0,0) 



(1,0) 



(0,1) 



(1,1) 



Figure 9: The possible perfect matchings in the case of conifold. Each perfect matching 
is associated with a lattice point of the toric diagram. 

3.3 Moduli space of vacua 

In this subsection, we study the moduli space A^D4 of the quiver quantum mechanics 



on the D4-D2-D0 state. The F-term conditions are now given by (3.2) and (3.3) 
together with 

dW 



dX n 







for 



X a 7^ Xp- 



(3.4) 



We particularly show that the moduli space is a subspace of the moduli space Aine 
of the parent D6-D2-D0 state obtained by replacing the D4-brane with a D6-brane. 
The subspace is characterized by the invariance under the action of a ?7(l)-subgroup 
of U(l) 2 x U(1)r. The [/(l)-subgroup depends on the divisor T> wrapped by the 
D4-brane. 

Perfect matchings 

To describe this, we first introduce so called "perfect matchings" of the bipartite graph 
Q y . Every edge in Q\ has its definite orientation from black to white vertex. A perfect 
matching m is defined as a collection of such oriented edges in Q y so that every vertex 
in Qq is attached to one and only one edge in m. All perfect matchings in the conifold 
case are shown in figure M We then define the "slope" (h x , h y ) of perfect matchings 
[T8l 152] . We first fix a reference perfect matching tuq and consider m a — itlq for every 
perfect matching m a , where the minus sign reverses the orientation of the edges, _ 
The edges in m a — m always form a closed oriented curve in Q v . When lifted to the 



16 To be more precise, we here regard m a and itiq as elements of a linear space over Z which is 
generated by the edges in Q v . For any edge e in Q v , we identify — e with the same edge of the 
opposite orientation. Then, for any two perfect matchings m a and m«, the element m a — mp can be 
identified with a closed oriented curve in Q v , or in T 2 . 
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Figure 10: The uplift of rri\ — m to the universal cover. The red arrows divide the 
universal cover into an infinite number of regions. The slopes are defined as the changes 
of the height function along the two periodic directions. This example has the slope 
(1,0). 



universal cover M 2 of T 2 such a closed curve divides 1R 2 into an infinite number of 



regions (figure 10). Then we define a height function on the universal cover, which 



takes a constant value in each region. We require that the height function changes by 
±1 when we move to an adjacent region crossing the curve. The sign of the height 
change depends on the orientation of the curve crossed. Then the slope (h x ,h y ) of 
m a is defined by the height changes in the two periodic directions of the quotient 



M 2 /Z 2 ~T 2 (figure [10). 

What is striking here is that h x , h y give a surjective map from the perfect matchings 
to the lattice points of the toric diagram E, where (h x , h y ) represents the relative 
positions of the lattice point [T71 [THl E2]- For example, (h x ,h y ) for all the perfect 
matchings in the conifold case are shown in figure [9] Note that the ambiguity of 
the choice of the reference matching mo is absorbed by shifting the origin. Since the 
lattice points of S are in one-to-one correspondence with toric divisors, we also have 
a surjection tp from perfect matchings to toric divisors. Although tp is not necessarily 
bijective, a divisor T>' associated with a corner of the toric diagram has a unique perfect 
matching m-D' such that ip(mx>>) = T>' if the brane tiling admits an isoradial embedding 
[oTll2"l] . Since our divisor T> is associated with a corner of the toric diagram, we have 
a unique perfect matching %. 
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Constraints on supersymmetric vacua 

We now claim that there are two constraints on the field configuration at supersymmet- 
ric vacua, which are derived from the F- and D-term conditions. The first constraint 



is 



J = 0, 



(3.5) 



which reduces (3.3) back to the original one dW /dX F = 0. Although (3.5) can easily 
be proven at a torus fixed point of the moduli space, it in fact holds at any point of 



the moduli space of vacua. We will prove this in subsection |3.7| The second constraint 
is that for any chiral field X a G Qi 

X a = if i/){X a ) G m v . (3.6) 

Here, recall that ip : Q — >■ Q w is a dual map. As we will see later, ip(Xp) is always 



included in the perfect matching %. Therefore the condition (3.6) is generically 



stronger than (3.2). However, we claim that (3.6) follows from (3.2) if combined with 



the 0-st ability for 9k < 0, 9* > and the other F-term conditions. We will prove this 
in subsection 13.81 



Now, the two different F-term conditions (3.2) and (3.3) from the D6-D2-D0 cases 



are simplified; the former is replaced with (3.6) while the latter reduces back to the 



original one. Moreover, the "anti-quark" J, which does not exist in the D6-D2-D0 case, 
vanishes on supersymmetric vacua. Hence, the moduli space «Md4 of the D4-D2-D0 



state is obtained just by imposing (3.6) on the moduli space A^D6 of the parent D6- 



D2-D0 state. Here the parent D6-D2-D0 state is obtained by replacing the D^-brane 
with a D6-brane. In other words, A^D4 is naturally regarded as a subspace of M.B6- 
This reflects the fact that our D4-brane extends only in a divisor of the Calabi-Yau 



three-fold while the D6-brane wraps on the whole three-fold. The constraint (3.6) is 
interpreted to mean that X a G il)~ l (mT>) describes a fluctuation of the fractional branes 
in a transverse direction to the D4-brane world-volume. As described in appendix [Cj 
we can verify this interpretation by considering a single DO-probe. 

From the mathematical point of view, our claim is the following. Let Q' be a new 
quiver such that Q' Q = Qo = Qo U {*} and Q[ — Q\ U {J} = Q± U {/, J}. Let us 
define T' to be the ideal of CQ' generated by dW/dl, dW/dJ and 8W/dX a for all 
A a GQi0 We then define 

A' := CQ'/T'. (3.7) 



"Recall that W = W + W &a 
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We also define an another algebra 

A:=CQ/F, (3.8) 

o ^ 

where T is the ideal of CQ generated by all the dW /dX a and X b such that ip(X b ) e 



rxij). Note that, since ip(Xp) G %, we have [Xp] = in A 18 Setting Xp = implies 
dW/dl = dW/dJ = and 8W/dX a = dW /dX a in CQ'. Therefore an A -module is 
naturally regarded as an A'-module. However, the converse is not always true. Now, 
our claim is that there is a one-to-one correspondence between ^-stable A'-modules 

o 

and ^-stable A -modules if we take 9k, 0* so that 9k < and 9* > 0. We will prove this 



in subsections 13.71 and 13.81 

M D4 in Mdg 

We here show that .M D4 as a subspace of A^D6 is characterized by its invariance under 
a [/(l)-subgroup of U{\) 2 x U(1)r ~ U(l) 3 . To see this, let us consider t e T such 
that 

[ 1 if ip(X a ) £ m v 

for a G 1. Since mp is a perfect matching of Q v , this t preserves all the F-term 
conditions for the parent D6-D2-D0 state. Furthermore, for general a, this t-action 
cannot be absorbed into gauge transformations. In fact, since each superpotential 



term includes one and only one chiral field involved in m-u, the t- action (3.9) does 
not preserve the superpotential. Such a [/(l)-action cannot be absorbed by gauge 
[7(T)IQo|-i_ Therefore, the t-actions (3.9) form a [/(l)-subgroup of U(l) 3 acting on 



Mm- We denote this by U{l) mv . 

Since all chiral fields involved in mx> vanish on M.^^ the moduli space A^D4 is 
invariant under the action of U(l) mv . Thus U(l) 3 reduces to U{\) 2 ~ U(l) 3 /U(l) mv 
on .Md4- The Witten index of the quiver quantum mechanics is then, up to sign, 
the number of L r (l) 2 -fixed points of A4d4- Since the perfect matching rrix> depends 
on the divisor T> wrapped by the D4-brane, so does the residual U(l) 2 . Note that all 



these properties of A^D4 essentially follow from the constraints (3.6) and (3.5). We 
will derive them in subsections 13.81 and 13.71 

Let us here mention that the relation between Aim and A^D6 is quite similar to 
that between the instanton and vortex moduli spaces. The moduli space of vortices 



18 Here, [Xp] is the equivalence class of Xp. 
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in d = 2, A/" = (4,4) theories was studied in [S3] and shown to be embedded in the 
moduli space of instantons|_] Moreover, it was pointed out that the vortex moduli 
space can be regarded as a U(l) invariant subspace of the instanton moduli space. It 
would be interesting to study this similarity further. 

3.4 Two-dimensional melting crystal 

We now consider the [7(l) 2 -fixed points of the moduli space .M D4- The inclusion map 
i : -Md4 ^ A^D6 implies that they are naturally regarded as some ?7(l) 3 -fixed points 
of -Md6- In other words, we are interested in a class of £/(l) 3 -fixed points of A^D6 
which are included in i(Air>i). Recall here that any £/(l) 3 -fixed points of A^D6 is 
expressed by a finite ideal tt of A*, where A* is the F-term equivalence class of paths 
starting at *. In particular, the elements of tt form the basis of the corresponding 
module. Now, what kind of finite ideal corresponds to a fixed point in i(A^D4)? The 



projection of «Md6 to \(M.m) is given by (3.6), which eliminates all the paths crossing 



the perfect matching mx>. Let us define A* as the set of F-term equivalence classes of 

o 

paths which do not cross the perfect matching mp. By definition, A* is a subset of 
A*. It is now clear that the ?7(l) 3 -fixed points in i(A^D4), or equivalently the U{1) 2 - 

o 

fixed points of -Md4> axe i n one-to-one correspondence with finite ideals of A*. From 

o 

the mathematical viewpoint, we find that [7(l) 2 -invariant ^-stable A -modules are in 

o 

one-to-one correspondence with finite ideals of A*. 

Next, recall that / : A* — > Ca„ is a bijection from A* to a three-dimensional crystal 
Ca,, where [ufcu/] G A* is mapped to an atom on the node k G Qq at the depth L 
Note that we here fix a reference node in Q and assume that the flavor-brane node * 
is attached to it; there is a quark / from * to the node. We denote by % the reference 

o 

node because it is an uplift of % G Qq. Let us now define C° := /(A*), which is a 

o 

subcrystal of Ca„- We call C° "reduced crystal." The finite ideals of A* are then in 
one-to-one correspondence with molten configurations of C° . The melting rule is the 
same as before; For a molten configuration p of C° , 

a bond from (5 G C° to a G p implies /3 G p. 

Note here that if one representative path of [vfco/] G A* includes some chiral field in 
^^mp) then so does any other representative path of [t^o/]. This implies that the 



19 See also [5H [SSJ [55] for the localization on the vortex moduli space of d = 2,J\f = (2, 2) theories. 
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Figure 11: Left: Two examples of zig-zag paths. In general, an edge in Q w is involved 
in two different zig-zag paths. There is a one-to-one correspondence between zig- 
zag paths and NS5-branes. For example, the red dotted line describes a NS5-brane 
associated with the red zig-zag path. Right: A zig-zag path V is decomposed into 
Zig(P) and Zag(P). 

reduced crystal C° is obtained from the original crystal Ca* by eliminating all the 
bonds associated with chiral fields involved in mp. 

Now, let us consider how the reduced crystal is embedded in the original three- 

o 

dimensional crystal. For any element [vkOJ e ] G A*, its representative path VkUr cannot 
contain any chiral field in if)~ l (mv)- However, u includes one such chiral field. Thus 

o 

we always have £ — 0, namely any element of A, is represented by a shortest path 
from *. This particularly implies that the reduced crystal C° lies in a two-dimensional 
plane at the depth zero. 

3.5 Shape of the crystal 

We have seen that the reduced crystal is always a two-dimensional crystal. We now 
discuss the precise shape of the crystal. Let us first introduce so-called "zig-zag paths." 
A zig-zag path is an oriented path in Q y which satisfies the following property: along 
a zig-zag path, we turn maximally right at a black vertex as well as maximally left 



at a white vertex (figure 11). This property is exactly the same as the property of 
cycles wrapped by NS5-branes in the brane tiling. In fact, there is a bijection from the 
NS5-branes to the zig-zag paths of Q v , which preserves the winding numbers. This 
also implies a one-to-one map between the zig-zag paths and the external legs of the 
toric web-diagram. The zig-zag paths of a dimer model with an isoradial embedding 
have some nice properties as described in appendix [Bj In particular, when lifted to the 
universal cover Q v , any two zig-zag paths share at most one edge. 
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Since a zig-zag path V is oriented, all the edges in V are also oriented. We then 
define Zig(P) as a collection of edges in V which are from white to black vertex. We 



also define Zag('P) := V — Zig('P) (the right picture of figure 11). Now, recall that 
our flavor brane is bounded by two NS5-branes at £±,£2- They are associated with 
two zig-zag paths in Q y , which we denote by V\ and V2 respectively. The winding 
numbers of Vi, V2 are determined by the directions of £±,£2- Then it was shown in |24j 
that 

Zig(Pi), Zag(P 2 ) C m v . (3.10) 

Note that, since V\ and V2 have different winding numbers, they always intersect with 
each other on Q y . At each intersection of V\ and V2, they share an edge ip(Y) for 



some chiral multiplet Y G Q\. Due to the condition (3.10), such a chiral multiplet Y 

always satisfies 

^{Y) G Zig(Pi) D Zag(P 2 ) C m v . (3.11) 



Note that Xp defined in subsection 3J2 is one such chiral multiplet. Therefore we 
always have ip(Xp) G mp. 

Let us now consider the uplifts of V\, V2 to the universal cover Q v , where we have 
a natural projection p v : Q v — > Q v . Recall that there is a reference node i in Qq ~ Q2 
which is attached to the D4-node *. The node i is also attached to a lift of the chiral 



multiplet X F , which we denote by X F (figure 12). Let ip : Q — >■ Q v be the dual map 
on the universal cover. There are two zig-zag paths of Q y which share the edge ip(Xp) 
at their intersection. When projected to Q v , they map to V\ and Vi- We therefore 
denote them by V\ and V2, respectively. The two zig-zag paths T > i,T > 2 divide the 
universal cover of T 2 into four regions. We denote one of them which includes % by 9L 



See figure 12 for example in the conifold case, where the blue face is % and the orange 



region including % is 9^. 

What we want to show is that the reduced crystal C° fills up the region 91 and 
V\ 1 V2 give the boundary of the crystal. Here, since C° lies at the depth zero, we 
identify atoms in C° with faces in Q y . Recall first that V\ and V2 share ip(Xp) at 
their intersection. The fact ip(Xp) G Zig('Pi) flZag(7 :, 2) then implies that Xp is always 



an out-going arrow from Dl and eliminated by the constraint (3.6 )r°] From this and 



20 Recall that our chiral multiplets have their definite orientations such that they encircle black 
(white) vertices in Q v clockwise (counter-clockwise). 
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Figure 12: The periodic dimer model Q y in the conifold case. The two zig-zag paths 
V\ and V<i divide the universal cover of T 2 into four regions. We denote by 9^ a region 
including the face %. We here assume that T> is associated with the upper-right corner 
of the toric diagram in figure [2j 



(3.10), it follows that all the out-going arrows from £H is eliminated by (3.6). This 



means that the reduced crystal C° is inside the region 91. 

Furthermore, if £ G Qo ~ Q^ is inside 91, then the reduced crystal C° always 
includes an atom placed on £. This can be shown as follows. Let us fix k, £ G Q and 
consider all the zig-zag paths of Q y which have k on its right side and £ on its left side. 
Such zig-zag paths are associated with some external legs of the toric web-diagram. 
We define C_ as a cone generated by such external legs. We also consider all the 
zig-zag paths on Q w which has k on its left side and £ on its right side, and define 
C + similarly. Then, it was shown in [21] that if the cone associated with our divisor 
T> is located in the clockwise direction of C+ and in the counter-clockwise direction 



of C_, as in figure 13, then there is always a path from k to £ which does not cross 
G pV )~ 1 ( m £>)r 1 | When we set k = i and £ to be inside 9\, then the cone associated with 
T> is always in the clockwise direction of C + and in the counter-clockwise direction of 
C_. Therefore, there is always a path from i to £ G Qo without crossing (p v ) _1 (mxi). 



21 Here (p v ) 1 (mx>) is the inverse image of m-p through p v : Q x 
matching of Q v . 



Q v , and itself is a perfect 
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Figure 13: If the cone associated with our divisor T> is in the clockwise side of C + and 
in the counter-clockwise side of C_, then there is a path from k to £ without crossing 
{p^)-\m v ). 



This means that C° includes an atom placed on £. Since £ is an arbitrary node in 9i, 
this implies that the reduced crystal C° fills up the region 91. 



3.6 Partition function 

We have shown that f/(l) 2 -fixed points of the moduli space A1d4 are in one-to-one 
correspondence with molten configurations of the reduced crystal C° lying in 9t. Then 
the partition function of the BPS D4-D2-D0 states is given by 

^D4-D2-D0 = ^(-l) dim ^) H x{\ (3.12) 

p feeQo 

where p runs over all possible molten configurations of C° and d t is the number of 
i-th atoms in p. We denote by .Mj-the moduli space of BPS states with charge {dk}, 
or equivalently the moduli space of vacua of the quiver quantum mechanics with ranks 
{dk} of the gauge groups. 

In general, the complex dimension of the moduli space dimc(.Mj*) is calculated as 

dim c (.A/f j*) =n 1 -n 2 - n 3 , (3.13) 

where n\ is the degrees of freedom of chiral fields which do not involved in the perfect 
matching m©, and n 2 is the number of non-trivial F-term conditions. In general rti 
can be written as 

wi = y] d s{X )dt{x), (3.14) 

X^m-p 



2<s 



where s(X) and t(X) is the starting and ending node of a chiral field X G Qi, respec- 



tively. The third term in (3.13) is given by 



n 3 =^(4) 2 , (3.15) 

fceQo 

which is the number of gauge degrees of freedom. Note that our gauge group has 

already been complexified to be ri/ceQo GL(dk), as explained in subsection 



2.3 



In the next section, we describe some examples and show that our construction 
perfectly reproduces known statistical models for D4-D2-D0 states on C 3 , (generalized) 
conifold, and C 2 /Ztv x C. We will also describe some examples of D4-D2-D0 crystal 
which have not been in the literature to the best of our knowledge. 

3.7 Proof of J = 

We here show that J = follows on a ^-stable module with 9k < for k G Qo- Since 
such a ^-stable module is a cyclic module generated by m G M*, it is sufficient to show 
that 

Jvl = (3.16) 



for any path v in Q. We assume that v is a path from i to j because otherwise (3.16) 
trivially holds. In general v can contain I and/or J. However, since I and J appear 
only as IJ in v, it is sufficient to consider v without including J, J. 

Hereafter, we consider the universal cover Q of Q and regard v as a path in Q. We 
respectively denote by i and j the starting and the ending nodes of v in Q, which are 
of course the uplifts of i and j. Now, for our purpose, it is sufficient to assume v is a 
shortest path from i to j. To see this, let us recall that our F-term equivalences are 

H; + IJ = °' w~ a = ° (for Xa * Xf) > (3 - 17) 

together with Xpl = and JXp = 0. If J = 0, these imply that v has a standard 
expression of the form [1] 

v = v co e °. (3.18) 

Here Vq is a shortest path from i to j, u is a loop which starts with Xp and surrounds 
a face in Q 2 , and £q is a non-negative integer. If J is non- vanishing, this expression is 
generally modified as 

v = v oo io + Wi{IJ)viU 1 ' 1 H h w n (IJ)v n uj en } (3.19) 
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Figure 14: A shortest path v s ry) from i to s(Y) for an out-going arrow Y from %\. 
Since Y crosses one of the red edges, there is a path v' such that Yv s (y) is F-term 
equivalent to v'Xp. 

where {v a } are shortest paths from i, and w a are paths to j which are not necessarily 
shortest. Here, by shortest path, we mean a path which cannot be decomposed further 
as (3.19). Since ul starts with Xpl = 0, the terms with £ a ^ do not contribute 
to Jvl. On the other hand, the terms with £ a = contribute Jw a IJv a I to Jvl. 
Therefore, if Jvl = is satisfied for all shortest paths v then it also holds for any 
other non-shortest paths v. 

Let us now show that a shortest path v which gets outside the region 9^ gives a 
vanishing contribution to Jvl. Recall that any out-going arrow Y from 9t is an element 
of Zig("Pi) U Zag(7 :> 2)- If v contains any such Y, then vl starts with Yv B ry\I where 
v s ry) is a shortest path from % to the starting node s(Y) of Y (figure 14). However, 



since there is always a path v' such that v'Xp is F-term equivalent to Yv s (y), we find 
Yv s (y)I = v'XpI = 0. Thus, if a shortest path v contains some out-going arrow from 
<K then Jvl = 0. 



The remaining task is to show (3.16) for shortest paths v which are inside the 



region 91. We first show the following lemma: 

Lemma 

A shortest path v from i to k e Q does not cross (p v ) _1 (mi)) if k is inside 9t. 



Proof. As mentioned in subsection 3.5, for any node k e Qq inside JH, there is a path 



v from i to k which does not cross (p v ) _1 (m£>). This path is, in fact, a shortest path 
from i to k. The reason for this is that any non-shortest path contain at least one u 
around some face in Q 2 , and therefore crosses (p v ) _1 (m,x)) at least once. □ 

This lemma particularly implies that a shortest path v from i does not contain any 
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Figure 15: The relative sign in (3.21 ) implies UiX F surrounds a black node while u 2 X F 
surrounds a white node in Q v . 



Xp if v is inside 9t. We now assume that v is such a shortest path. We rewrite Jvl as 

= tr (mv) - tr (u 2 v) , (3.20) 



Jvl = tr(IJv) = — tr ( ^„ /' 



dX F 
where U\ and w 2 are paths in Q such that 

dXp" ~ 



Mi + U 2 . 



(3.21) 



The relative sign implies that U\Xp surrounds a black node in Qq while u 2 Xp sur- 



rounds a white node (figure 15). Note that U\V and u 2 v form loops in Q, which we 



denote by L\ and L 2 respectively. In terms of these, we can rewrite (3.20) as 

Jul = tr(£i) -tr(L 2 )- 



(3.22) 



In the rest of this subsection, we will show that this vanishes for any shortest path v 
inside 9t. 

Rhombus tiling and F-term equivalence 

Note that the two loops L\, L 2 have the same winding number (n x , n y ). Since v is inside 
the region 01, {n x ,n y ) always lies in a convex cone Cx> in K 2 . The cone is naturally 



identified with J defined in subsection |3.2| and therefore depends on the choice of V. 
We can generally write the winding number as 



\Tlx) fly) 



r{n Xl n y ), 



(3.23) 



where r is a positive integer and n x ,n y are mutually prime integers. Since r depends 
on the shortest path v, we sometimes write it as r v . 

Recall here that there is a one-to-one correspondence between zig-zag paths and 
the external legs of the toric web-diagram. In particular, if there is a zig-zag path 
whose winding number is (n' x ,n' y ) then there is an external leg spanned by the vector 
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Figure 16: The situation in the case of r = 1. The blue dashed rectangle expresses 
the fundamental domain for the torus IR 2 /Z 2 . The black solid edges together with 
the while and black nodes form the periodic dimer model Q w . The union of the gray 
rhombi is denoted by Hp. When we regard L\, L 2 as paths in IR 2 , they are continuously 
deformed into each other without crossing Hp. 



(n' x , n'y). Since 5 is a cone spanned by two adjacent external legs £%, £ 2 , the only zig-zag 
paths whose winding number lies in Cx> are V\ and V 2 . This fact will be important in 
the proof of tr(Li) — tr(L 2 ) = 0. 

The case r = 1 

We first consider the case r = 1. We will show that L% and L 2 are F-term equivalent 



to each other, which implies (3.22) vanishes. Recall that our brane tiling admits an 
isoradial embedding, and therefore we have a rhombus tiling on T 2 . We denote by Hp 
a rhombus including the chiral multiplet Xp. What is important here is that, when 
we regard Li,L2 as cycles of T 2 rather than those of Q, they can be continuously 
deformed into each other without crossing Hp. To see this, we move to the universal 
cover IR 2 . The uplifts of Li,L 2 are paths with infinite length, which we denote by 
Li and L 2 respectively. Note that each of L 1; L 2 has an infinite number of connected 
components, which are just copies of one connected line with infinite length. Since v 
is a shortest path in Q, they have no self-intersection. The uplift of the rhombus Hp 
is an infinite number of rhombi. We denote by Hp the union of all such rhombi in IR 2 . 
We now see that L\ can be continuously deformed into L 2 without crossing Hp (figure 
16). This always follows when r = 1. Furthermore, we can perform this deformation 
periodically so that L\ remains the uplift of some cycle of T 2 . This implies that, on 
T 2 , Li can be deformed into L 2 without crossing H F . 
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Figure 17: The collection of rhombi passed by Li,L 2 (white rhombi) divides the 
universal cover into an infinite number of connected regions. The union of the regions 
which do not contain any Xp is denoted by x. In this picture, x is the union of the 
orange rhombi. 



Note that this does not immediately mean L\ and L 2 are F-term equivalent. When 
we regard them as paths in Q, rather than those in T 2 , we can deform them only by 
F-term equivalence. In fact, in a general dimer model, there could exist two homotopic 
paths which are not F-term equivalent to each other. However, there is a nice theorem 
that any two homotopic paths are F-term equivalent if the dimer model admits an 
isoradial embedding [57]. Below, we will generalize this and show that L x and L 2 are 
in fact F-term equivalent. 

In the universal cover IR 2 , there are special rhombi which are passed by L\ or L 2 . 
The collection of such rhombi divides IR 2 into an infinite number of connected regions. 



We denote by r the union of the regions which do not contain any Xp (figure 17). Note 



that the boundary of r has two kinds of vertices; quiver vertices and dimer vertices. Let 
us now consider to reduce the region x via the F-term equivalence. We deform L\ by 
using the F-term equivalence at the boundary of r, which reduces the area of x (figure 
18). Reducing x is possible only when the boundary of x has a dimer vertex which is 
not attached to any rhombus edge lying inside x [57] . We perform this deformation of 
L\ 'periodically so that it is an uplift of a deformation of L\. Then this does not change 
the winding number of L\. We repeat this deformation and stop if we cannot reduce 
x further. There are, in fact, only three possibilities to stop. The first possibility is 
that x is empty when we stop. In this case L\ coincides with L 2 after the operations, 
and therefore we find that L\ is F-term equivalent to L 2 . The second possibility is 
that, when we stop, r is not empty but can be written as the disjoint union of compact 
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Figure 18: By using the F-term equivalence at the boundary of the region t, we can 
reduce the area of r. We here used the F-term condition associated with the blue edge. 
Reducing r is possible only when we have a dimer vertex at the boundary of r which 
is not attached to any rhombus edge lying inside r. 
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Figure 19: Left: The second possibility to stop. When we stop, the region r is not 
empty but written as the disjoint union of compact regions. Right: The third pos- 
sibility to stop. When we stop, r is not empty and cannot be written as the disjoint 
union of compact regions. 



regions. In this case, we have the situation illustrated in the left picture of figure 19 



However, this is inconsistent with the rhombus tiling of IR 2 , which follows from Lemma 
5.3.1 of [57]. Thus, the second possibility is forbidden. The third possibility is that, 
when we stop, r is not empty and cannot be written as the disjoint union of compact 



regions. In this case, we have the situation as in the right picture of figure 19 This 
does not contradict with the lemma in |57| . but still turns out to contradict with the 
rhombus tiling of IR 2 . To see this, note first that any dimer vertex at the boundary of 
r is now attached to at least one rhombus edge lying inside r, because otherwise we 
can further reduce r. It also follows that any quiver vertex at the boundary of r is 
attached to at least one such rhombus edge. Therefore the simplest case is that, at the 
boundary of r, each (quiver or dimer) vertex is attached to only one edge lying inside 



tr. In this case, there is always a zig-zag path V included in the region r (figure 20). 
When projected to T 2 , V is mapped to a zig-zag path of Q whose winding number is 
(n x ,n y ). Since (n x ,n y ) lies in the convex cone Cx>, such a zig-zag path should be V\ 
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Figure 20: The simplest case of the third possibility. There is a zig-zag path inside x 
(blue path) which has the same winding number as L\. In this situation, we cannot 
deform L\ to reduce x because every white dimer vertex at the boundary is attached 
to some rhombus edge lying inside r. 



or V%. Recall that both V\ and V 2 involve Xp. Therefore V also involves Xp. This, 
however, contradicts with the fact that the region x does not contain any Xp. Hence, 
the simplest case does not occur. Let us then consider a non-simplest case. Namely, 
suppose that the boundary of t includes a vertex attached to more than one rhombus 
edge lying inside x. In this case, L\ has to intersect with itself in order to be periodic 
in IR 2 . The periodicity of L\ is necessary because it is the uplift of L\. However, the 
self-intersection of L\ contradicts with the fact that v is a shortest path in Q. This 
implies that not only the second possibility but also the third possibility is forbidden. 



Hence, L\ is always F-term equivalent to L 2 , which implies (3.22) vanishes. 



The case r > 2 

Next, we turn to the cases r > 2. In this case, Jvl = follows by induction on 
r. Since r depends on v, we here write it as r v . We assume Jv'I = for any v' with 
r v i < r v , and prove Jvl = 0. 

Note first that, in the case r v > 2, it is not possible to deform L\ into L2 without 
crossing Hp. A typical example with r v = 3 is shown in figure 21 However, as we 



will prove below, there are always one-cycles L\^L ! 2 in Q which can be continuously 
deformed into each other without crossing Hp and also satisfy 

tr(Li) - tr(L 2 ) = tr^) - tr(L' 2 ) - (r v - 1) tr(vIJ). (3.24) 

Once such L[, L' 2 exist, we can show 

r v (JvI) = tr(Li) - tr(L' 2 ). (3.25) 

Since L[ can be continuously deformed into L' 2 without crossing Hp, the same argu- 



ment as in the case of r v = 1 implies that (3.25) vanishes. Since r v > 0, we then find 
Jvl = 0. 
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Figure 21: When r > 2, it is not possible to deform L% into L 2 without crossing Hp. 
The picture shows the universal cover in the case of r = 3. The blue dashed rectangle 
expresses the fundamental domain for the torus IR 2 /Z 2 . 

The remaining task is to show the existence of such a nice set (L' t , Z/ 2 ). Let us first 
consider an example in figure 21, where Li,L 2 are lifted to Li,L 2 , respectively. We 

Since they are periodic in the universal 
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now define paths L[ and L' 2 as in figure 
cover, we can regard them as the uplift of some one-cycles of T 2 . We denote such 
one-cycles by L[ and L 2 , respectively. By definition, L[ and L 2 can be continuously 
deformed into each other without crossing Hp. Let us now consider the difference 
between L\ and L' x . In fact, the only difference between them is that, only once, they 



pass the opposite side of H F . In terms of u\, u 2 defined in (3.21 ), this can be expressed 



as 



tr(Li) — tr(L / 1 ) = tr {v'u2) — tr(t>'wi), 



(3.26) 



for some path v' from i to j (figure 23 ). The relation (3.21 ) and the F-term equivalence 



(3.3) now imply that this can be rewritten as 

tr(L 1 )-tr(L / 1 ) = -tr(v , U). 



(3.27) 



Note here that v' and v are different paths, but v'U and vIJ have the same winding 
number r v (n x ,n y ). In fact, the only difference between v and v' is again that they pass 
the opposite side of Hp once. Then, the same argument as above implies that there 



are paths v i , v 2 (figure 23 ) from i to j such that 



tx{v'IJ) - tx(vIJ) = tr(t>x Uv 2 IJ). 



(3.28) 



Here the right hand side can be written as (Jv\I)(Jv 2 I). Note that, since v\ and v 2 
are shorter than v, we always have r Vl ,r V2 < r v . Then, due to our assumption, both 
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Figure 22: There are always L\ and L' 2 which can be deformed into each other without 



crossing Hp and also satisfy (3.24). The picture shows the uplifts of (Li,L 2 ) for 



(Lx, L 2 ) in figure 21 The blue dashed rectangle expresses the fundamental domain for 
the torus. 



Jv\I and JV2I vanish. This implies that we can rewrite (3.27) as 



tr(Li) - tr(L;) = -tr(vIJ). 



(3.29) 



The minus sign in the right-hand side is important, which comes from the fact that 
L\ always surrounds the black-node side of Hp. Let us now turn to L 2 and L' 2 . The 
similar argument as above tells us that we always have 



tr(L 2 ) - tr(l/ 2 ) = +tr(vIJ). 



(3.30) 



The different sign in front of tr(vIJ) is due to the fact that L 2 always surrounds the 



white-node side of Hp. Combining (3.29) and (3.30), we obtain 



tr(Lx) - tr(L 2 ) = tr(L' x ) - tr(L 2 ) - 2tr(vIJ). 



(3.31) 



This is a special example of (3.24) with r v = 3. 



By generalizing the above argument, we can show that the existence of the desired 
(L[, L' 2 ) for any shortest path v which lies inside 9\ with r v > 2. The crucial point is 
that the winding number (n x , n y ) of L%, L2 for such v always lies in the convex cone Cv- 
In general, we can deform L 1; L 2 into the following forms by the F-term equivalence 
without crossing H F : 



L\ ~ V (1) U S1 ■ ■■V {rv ^ 1) U Srv _ 1 V( rv) U 1 



V(1)U S1 ■ ■ ■ V(r v -l)U arv _ x V {rv )U 2 , (3.32) 
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Figure 23: In the example of figure 21 there is a path v' from i to j such that 



tr(Li) = tr(t>'u 2 ) and tr(L' 1 ) = tr(v'ui). In the picture, the red path is v' while the 



blue paths are U\ and w 2 . The green V\ and t> 2 are used in (3.28). 



where v^) is a shortest path with r v = 1, and Sk G {1, 2}. Then our desired L[,L' 2 
are defined by 



L 'l = W (1) U 1 • " " V (r v -l)UlV(r v )Ul, L' 2 = V(l)U2 " " " ^-l)^ (r„) W 2- 



(3.33) 



The fact that (n x ,n y ) lies in the convex cone Cx> guarantees that L[ can be deformed 
into L' 2 without crossing Hp. The difference between tr^i) and tr(L' 1 ) is generally 
expressed as 



tr(Li) - tr(Li) = -nitr(v/«7), 



(3.34) 



where n\ = J2k=i ( s fe — 1) > 0. Here, we used our assumption that Jv'I = for any 
shortest path v' with r v > < r v . Similarly, the difference between tr(L 2 ) and tr(L 2 ) is 
written as 



tr(L 2 ) - tr(L 2 ) = +n 2 tr(vIJ), 
where n 2 = Y7k=i (^ — s k) > 0. Then we obtain 

tr(L0 - tr(L 2 ) = tr(L;) - tr(L 2 ) - (r„ - l)tr{vIJ), 



(3.35) 



(3.36) 



where we used the fact that n\ + n 2 = r v — 1. Thus, we have shown the existence of 
the desired (L' 1; Z/ 2 ), which completes our proof of J = 0. 



3.8 Proof of m p = 



We here prove (3.6) on a ^-stable module with 9k < for k £ Q\. Since such a ^-stable 



module is a cyclic module generated by an element m G M#, it is sufficient to show 
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that ifYEQi satisfies (ip o p)(Y) G m-p then Yvl = for any path v from i G Q 
According to the argument in the previous subsection, we here set J = 0. 

Let s(Y) be the starting node of Y. We assume v is a path from % to s(Y) because 
otherwise Yvl = trivially holds. Moreover, it is sufficient to consider only the 
shortest path v = v s ry) from % to s(Y). To see this, recall that v has the standard 



form (3.18) for the shortest path vq = v s (y) and a loop uj starting with Xp. As shown 



in [3], this can also be expressed as 

v = u e v s (Y)i (3.37) 

where a) is a loop which starts at s(Y) and surrounds a face in Q2. In particular, we 
can take u to start with Y. Then Yvl = immediately follows if Yv s iy\I = 0. Hence, 
we only need to consider the case v = v s (y) in our proof. 

The remaining task is to prove Yv s (y)I = for Y G Q\ satisfying (ip op){Y) G m-p. 
Let us first consider FeQi which is an out-going arrow from 9t. Recall that such Y is 
always an element of Zig("Pi) U Zag(7 :> 2), and therefore satisfies (ip o p)(Y) G %. The 
path Yv s (y)I for such Y has already been considered in the previous subsection, and 
shown to vanish. This also implies that any Y G Q\ outside 9i leads to Yv s (y)I = 0. 

Now, let us consider Y G Q\ which is inside %K and satisfies (ip ° p){Y) G mp. 
We here denote by t(Y) the ending node of Y. Since t(Y) is inside 1H, our lemma in 



subsection 3.7 implies that the shortest path v t {y) does not cross (p v ) 1 (m v ). Now, 



let u>t(Y) be a loop around a face in Q which starts at t(Y) and ends with Y. We then 



define a path uy so that cot(Y) = Yuy (figure 24). Here, the path Uyv t (y) is a path 
from % to s(y) without crossing (p v ) _1 (mx)), which implies that UyV t (y) is a shortest 
path to s(Y). We can then write u«(y) = Wy^t(y) up to the F-term equivalence, and 
therefore 

Ff s (y)J = Yu Y V t (Y)I = U t {Y)V t {Y)I ■ (3.38) 

As mentioned already, Wt(y) v t(y) is F-term equivalent to v t (y)0J where w is a loop 
starting with Xp. Since v t (y)0jl starts with Xpl, it vanishes due to the constraint 
(3.2). Thus we find Yv s iy\I = 0. Combining the argument in the previous paragraph, 



we have shown Yv s (y)I = for Y which is inside 9^ and satisfies (ip o p)(Y) G m-p- 
This completes our proof. 



22 Recall here that ifj '■ Q —> Q V is the dual map. 
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Figure 24: There is always a path from i to t(Y) without crossing (p v ) _1 (mp) if Y is 
inside 91. Then Uyv^y) is a shortest path from i to s(V). Here the red edges are the 
ones in (p s/ )~ 1 (rri'p). 

4 Examples 

In this section, we consider some examples of our crystal melting model for D4-D2-D0 
states on toric Calabi-Yau three-folds. 



4.1 C 3 

The first example is C 3 , where we can put a single D6-brane on the whole C 3 or a 
D4-brane on a divisor C 2 . In either case, we have no D2-brane charge because C 3 has 
no compact two-cycle. 

Let us first consider the D6-D0 bound states. The toric diagram and the brane 



tiling of C 3 are shown in figure 25 In the quiver quantum mechanics on the DO-branes, 
we have one gauge group and three chiral multiplets B a . The rank of the gauge group 
is equal to the DO-charge. The quiver diagram of the theory is shown in the right 



picture of figure 25 From the brane tiling, the superpotential is read off as 



W = ti(B 1 [B 2 ,B 3 ]) 
which implies the F-term condition 

[B a ,B b ] = 0. 



(4.1) 



(4.2) 



This implies that any element of A* is now expressed as a monomial of the form 
(Bi) ni (B 2 ) n2 (B 3 ) n3 I for rii G N. The periodic dimer model Q y is depicted as in 
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Figure 25: Left: The toric diagram of C 3 . The lattice points are at (0,0), (1,0) and 
(0, 1). Middle: The brane tiling for DO-branes on C 3 . Right: The quiver diagram for 
D6-D0 states on C 3 . 



the left picture of figure 26 Note that this dimer model clearly admits an isoradial 
embeddingr 3 ! By putting atoms on the faces of Q v , we obtain a cubic crystal as in 



the right picture of figure [26j where a single cube is an atom associated with the unit 
DO-charge. Each face of the cube expresses a "bond" associated with a chiral multiplet 
B a . The melting rule of the crystal implies that the partition function (2.9) is now 
written as a sum over three-dimensional Young diagrams p [21 E] : 



p n =l V y ' 

where q is the Boltzmann weight for DO-charge. 

Now, instead of the D6-brane, let us put a D4-brane on a holomorphic divisor C 2 in 
C 3 . We consider BPS DO-branes bound to the D4-brane. Recall that this example has 
already been mentioned in the introduction. The perfect matchings of Q w are shown 



in figure [27} We have three perfect matchings, corresponding to the three corners of 
the toric diagram of C 3 . Let us assume, without loss of generality, that the divisor 
T> wrapped by the D4-brane is associated with the lattice point (0, 1) of the diagram. 
Then, the boundary NS5-branes of the facet 5 are the second and third NS5-branes in 



figure 25 Since Bj, is located at their intersection, we identify Xp = B% in this case. 
The superpotential induced by the D4-brane is then written as 



W &a 



JB«I. 



(4.4) 



The perfect matching mx> only involves B%, and therefore the condition (3.6) just 
implies B% = on supersymmetric vacua. Setting S3 = in fact leads to the well- 



2:ic 



See appendix Ib] for a general criterion. 
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Figure 26: Left: The universal cover Q y of the dimer model Q y for C 3 . Each dia- 
mond region is a fundamental domain of IR 2 /Z 2 ~ T 2 . The two zig-zag paths which 
are attached to i and share B 3 at their intersection are also shown. Right: The 
corresponding three-dimensional crystal is obtained by putting atoms on the faces of 



known D4-D0 quiver on C 2 (figure 28). The only non-trivial F-flatness condition is 
now 



[B 1 ,B 2 ] + IJ = 0, 



(4.5) 



which is equivalent to the ADHM constraint. Thus, we have reproduced the well- 
known result of the D4-D0 states on C 2 . Note that we already know that J = also 
holds on supersymmetric vacua. 

The [7(l) 2 -fixed points of the moduli space are in one-to-one correspondence with 
the molten configurations of a two-dimensional crystal in figure [28j where a single 
square atom is associated with the unit DO-charge. The four sides of each square 
express four "bonds" attached to it. Note that the crystal is a two-dimensional "slope 



face" of the three-dimensional crystal in figure [26J Moreover, the boundary of the 
two-dimensional crystal is given by the two zig-zag paths shown in the left picture of 



figure 26 The melting rule of the crystal implies that every molten configuration is 



expressed as a two-dimensional Young diagram. The partition function (3.12) is then 
written as a sum over two-dimensional Young diagrams p: 



oo 1 



We here implicitly assume that all the fixed points are bosonic. We can show this by 
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Figure 27: Leftmost: The dimer model Q y on T 2 for C 3 . Others: Three perfect 
matchings corresponding to the three lattice points of the toric diagram. In C 3 case, 
the perfect matchings are in one-to-one correspondence with the toric divisors of C 3 . 
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*cwo 



IB, 




Figure 28: Left: The quiver diagram for D4-D0 states when T> is associated with the 
lattice point (0, 1) of the toric diagram. Right: The two-dimensional crystal for the 
D4-D0 states on C 2 . 



calculating the dimension of the moduli space. The moduli space is parameterized by 
/, J, B\ and B 2 which have 2|p| + 2|p| 2 complex parameters. They are subject to the 
F-term conditions [Bi, B 2 ] +IJ = which reduces |p| 2 of them. Furthermore, dividing 
out the gauge degrees of freedom reduces further |p| 2 complex parameters. Thus, the 
complex dimension of the moduli space is 



dim c (M 



IpI, 



2|P|, 



(4.7) 



which implies that all the fixed points of the moduli space are bosonic. Note that (4.7) 
matches the dimension of |p|-instanton moduli space. 



4.2 Conifold 

The second example is the conifold. We have already seen the brane tiling, D6-D2-D0 
quiver, and D6-D2-D0 crystal in section |2j The partition function associated with 
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Figure 29: Left: The quiver diagram for D4-D2-D0 states on the conifold. Right: 
The reduced crystal is a two-dimensional slope face of the three-dimensional pyramid 
crystal in figure [6j 



the three-dimensional melting crystal in figure p\ was considered in [3] in the study 
of the non-commutative Donaldson-Thomas invariants. In this subsection, we instead 
consider D4-D2-D0 states on the conifold. Note that the dimer model Q y for the 
conifold in figure [5] admits an isoradial embedding. 

We first put a flavor D4-brane on a toric divisor T> associated with a corner of the 
toric diagram in figure [2] Since all the toric divisors are isomorphic in the singular 
conifold limit, we assume without loss of generality that T> is associated with the 
lattice point (1,1) of the toric diagram. The boundary NS5-branes are then the first 
and second ones in the right picture of figure [3j Since B 2 comes from the massless 
string at their intersection, Xp is identified with B 2 in this case. The superpotential 
induced by the D4-brane is written as 



^flavor = JB 2 I. 



(4- 



The total superpotential is given by W + Wfl avo r with 



W = ti(A 1 B 2 A 2 B 1 ) - ti(A 1 B 1 A 2 B 2 



(4.9) 



The perfect matching mp turns out to be 7713 in figure 9 According to (3.6), placing 



a D4-brane on T> eliminates all the chiral fields involved in mp from the massless 



spectrum, leading to a quiver diagram shown in the left picture of figure 29 The 
non-trivial F-term condition is now 



A 2 B 1 A l - A X B X A 2 + IJ = 0. 



(4.10) 



The two-dimensional crystal for the D4-D2-D0 states then extends in the region 91 in 



figure 12 As shown in the right picture of figure 29 , this crystal is a two-dimensional 
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Figure 30: Left: The two-dimensional crystal for D4-D2-D0 states on the conifold, 
which infinitely extends in the upper-right region. Three chiral fields Ai,A 2 ,B 2 cor- 
respond to three types of "bonds" connecting triangle atoms. Right: An example of 
a molten crystal. 

"slope face" of the three-dimensional pyramid crystal in figure |6j Note here that we 
already know that J = also follows on the supersymmetric vacua. 

Let us express the two-dimensional crystal in a different way, keeping the crystalline 
structure. We express each atom as a triangle so that its three sides are associated 
with the three "bonds" attached to the atom. Then the two-dimensional crystal in 



figure [29] is re-expressed as in the left of figure 30 This two-dimensional crystal is 
precisely equivalent to that of the triangular partition model proposed in [T3]. The 
melting rule of the crystal implies that 

• A blue triangle can be removed if and only if its left and lower edges are not 
attached to other atoms. 



• A red triangle can be removed if and only if its slope edge is not attached to 
other atoms. 

These are exactly the same rules as proposed in [13J. 

The partition function of the BPS index is written as a sum over the molten 
configurations of the crystal: 



Crystal == ^(-l)^^.^V, 



(4.11) 



where p runs over all the molten configurations, and a and b are the numbers of blue 
and red triangular atoms in p. The sign factor is determined by the dimension of 
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the moduli space. We now have five non-vanishing chiral fields J, J, Ai, A 2 , B\ with 
3ab + a + b complex degrees of freedom in total. The non-trivial F-term condition 



(4.10) reduces ab of them. Dividing out the gauge degrees of freedom further reduces 



a 2 + b 2 parameters. Thus, the dimension of the moduli space is written as 

dim c {M a , b ) = a + b-{a-b) 2 . (4.12) 

This implies that 

(_ 1 )dfan C (M , b ) = 1} (4 13 ) 



and therefore the fixed points of the moduli space are all bosonic ._ When we change 
the variables as 

q = xy, Q = -x, (4-14) 



the partition function (4.11) can be rewritten as [13] 

oo 1 oo 

n i3^ n (!-«)• ( 4 - is ) 



-^"crystal 



q 

n=l m=0 



By identifying q and Q with the Boltzmann weights for DO and D2-branes respectively, 
this corresponds to the correct BPS D4-D2-D0 partition function in the singular limit 
of the conifold 



4.3 Suspended pinch point 

The third example is the suspended pinch point, whose toric diagram is shown in the 



left picture of figure [31] In this case, we have four toric divisors associated with the 
corners of the toric diagram. We first consider a D4-brane on the divisor for pi, and 
next consider that for p 2 . The former reproduces the oblique partition model proposed 
in [H] while the latter gives a new statistical model. It will be shown via the wall- 
crossing formula that they give the correct partition function of the BPS indices for 
the corresponding D4-D2-D0 states. 

D2-D0 Quiver and bipartite graph 

We first identify the quiver diagram for the D2-D0 states on the suspended pinch point. 



The relevant brane tiling is shown in the middle of figure 31 The quiver diagram of the 



24 Note here (4.12| is always non-negative. In fact, our melting crystal implies b > — — — ' ""'""" ' ' 



which is equivalent to a + b > (a — b) 
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Figure 31: Left:The toric diagram of the suspended pinch point. Right: The brane 
tiling for the D2-D0 states on the suspended pinch point singularity. 

theory on the D-branes is shown in the right picture. We have seven chiral multiplets 
with the following superpotential: 



W, 



= tr ($(YY - XX)) + tr [zZXx\ - tr (zZYY\ . 



(4.16) 



The F-term conditions are then written as 



YY = XX, ZXX = YYZ, XXZ = ZYY, 
$X = XZZ, X$ = ZZX, $Y = YZZ, Y$ = ZZY. (4.17) 

If we introduce a flavor D6-brane, then we have a D6-node attached to one of 
the nodes in the quiver diagram without any additional superpotential. The U(l) 3 - 
fixed points of the moduli space are expressed as molten configurations of the three- 
dimensional crystal studied in [5 J (In particular, see figure 5). 

In this subsection, we instead introduce a flavor D4-brane on a toric divisor. The 



dimer model Q and its perfect matchings are shown in figure 32 Note that Q admits 
an isoradial embedding. 



Oblique partition model 

We first put a D4-brane on a divisor associated with the lattice point p\ of the toric 
diagram. The boundary NS5-branes for the divisor are the first and the second NS5- 



branes in figure 31 Since the chiral multiplet Z comes from the massless string at the 
intersection of the boundary NS5-branes, it is Z that should be identified with Xp in 
this case. This implies that our D4-node * is now attached to the blue and red nodes 
in the quiver diagram. The additional superpotential induced by the D4-brane is now 
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Figure 32: The dimer model Q y and perfect matchings for the suspended pinch point 
singularity. The coordinate below each perfect matching expresses the corresponding 
lattice point of the toric diagram. 



written as 



Wf 



flavor 



JZI. 



(4.18) 



As shown in figure 32, the perfect matching associated with pi is m 2 , which includes 



Z and $. The constraint (3.6) on supersymmetric vacua is now written as 

z = 0, $ = 0. 



(4.19) 



This means that Z and $ become massive with vanishing vev's in the infrared. The 
quiver diagram from which the two massive fields are eliminated is shown in the left 



picture of figure [33} The non-trivial F-term conditions are now 
YY -XX = 0, ZXX - YYZ + IJ = 0. 



(4.20) 



Note here that, from the general argument in the previous section, we already know 
that J = also follows on the supersymmetric vacua. We here explicitly write J for 
comparison with the previous C 3 example. 
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Figure 33: Left: The quiver diagram for the D4-D2-D0 states when the D4-brane is 
on a divisor associated with p\ in the toric diagram. The D4-node is denoted by *. 
Right: The two-dimensional crystal for the D4-D2-D0 states, which reproduces the 
so-called "oblique partition model" proposed in 



The two-dimensional crystal for the D4-D2-D0 states is expressed as in the right 



picture of figure [33j where each atom is assigned the same color as the corresponding 
quiver node. The (three or four) sides of each atom are associated with "bonds" 
attached to it. This crystal is the same as that in the "oblique partition model" 
proposed in [T4]. Our melting rule of the crystal implies the followings: 

• A blue triangle can be removed if and only if its left and lower edges are not 
attached to other atoms. 

• A red triangle can be removed if and only if its slope edge is not attached to 
other atoms. 

• A parallelogram can be removed if and only if its left and lower edges are not 
attached to other atoms. 



These are exactly the same rules as proposed in [14] . Note that, as mentioned in the 
previous section, this two-dimensional crystal is a "slope face" of the three-dimensional 
crystal associated with the parent D6-D2-D0 states. 

The partition function of the BPS index is evaluated by counting molten configu- 
rations p of the two-dimensional crystal: 



^crystal = ^(-l)^^'- VyV, 



(4.21) 
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where a, b, c are the numbers of blue triangles, red triangles and white parallelograms 
in p, respectively. The variables x, y, z are the Boltzmann weights for the atoms. The 
sign factor is determined by the dimension of the moduli space Ai a ,b- We now have 
seven non-vanishing chiral fields with (a+b+ab+2bc+2ca) complex degrees of freedom 



in total. The non-trivial F-term conditions (4.20) reduce (c 2 + ab) of them. The gauge 



degrees of freedom reduce further (a 2 + b 2 + c 2 ) parameters. Thus the dimension of 
the moduli space is given by 

dim c {M aAc ) = a + b + 2bc + 2ca- 2c 2 - a 2 - b 2 . (4.22) 

Then the sign factor is given by 

r^\dim c {M a ,b,c) = i (4.23) 

In terms of different variables 

q = xyz, Q 1 = -x, Q 2 = z, (4.24) 

the partition function is written as 

oo _ oo 

Crystal = J] 7"^ II i 1 ~ ^M 1 ~ *"<*&*)' C 4 " 25 ) 

n=l m=0 

It was shown in [T4J via wall-crossing formula that this is the correct BPS partition 
function of the D4-D2-D0 states on the suspended pinch point singularity, where q is 
the Boltzmann weight for DO-charge and Qi,Q2 are those for D2-charges associated 
with two blowup two-cycles. 

Another statistical model 



We now put a D4-brane on the divisor associated with p 2 in figure 31 The boundary 



NS5-branes are then the third and fourth ones in figure 31 Since Y comes from the 
string at the intersection of the boundary NS5-branes, Xp is now identified with Y. 
This implies that the D4-node * is attached to the blue and white nodes in the quiver 
diagram. The superpotential induced by the D4-brane is now written as 

W &avoI = JYI. (4.26) 



The perfect matching associated with the lattice point p 2 is mo in figure 32 The 



constraint (3.6) then implies that 

X = 0, Y = 0, (4.27) 
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Figure 34: Left: The quiver diagram of the theory on D4-D2-D0 states with a D4- 
brane on a divisor associated with P2 in the toric diagram. Right: The corresponding 
two-dimensional crystal. 



in the infrared. After eliminating X and Y, the quiver diagram can be depicted as in 



the left picture of figure 34 The non-trivial F-term conditions are now 
$X - XZZ = 0, Y$- ZZY + IJ = Q. 



(4.28) 



Although J = also follows in the infrared, we here explicitly write it for comparison 
with the other examples. 



The two-dimensional crystal can be depicted as in the right picture of figure [34 
Note that this crystal is an another "slope face" of the same three-dimensional crystal 
for the parent D6-D2-D0 states. In fact, choosing a different divisor T> of the same 
Calabi-Yau leads to a different slope face of the same three-dimensional crystal. The 
melting rule of the crystal now implies 

• A blue triangle can be removed if and only if its left and lower edges are not 
attached to other atoms. 

• A red triangle can be removed if and only if its slope edge is not attached to 
other atoms. 

• A white rectangle can be removed if and only if its left and lower edges are not 
attached to other atoms. 

This two-dimensional melting crystal model is, to the best of our knowledge, a new 
example of statistical model for D4-D2-D0 states. We therefore examine it in detail. 
The partition function of the BPS index of the D4-D2-D0 states is given by 



2 crystal = ^(-lf^^M^y^, 



(4.29) 
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Figure 35: Left: A molten configuration of the crystal. Right: The left configuration 
is split into towers of atoms. Each tower has several squares composed of a blue 
triangle, a red triangle and a white rectangle. Some of the towers have an additional 
blue triangle with or without a red triangle on its top. 

where a, b, c are the numbers of blue triangles, red triangles and white rectangles in 
p, respectively. The variables x,y and z are Boltzmann weights for the three types of 
atom. The sign factor is determined by the dimension of the moduli space Ai' a bc . We 
here have seven massless chiral fields with (a + c + 2ab + bc + ca + c 2 ) complex degrees 
of freedom in total. The F-term conditions ( |4.28[ ) reduce (bc + ca) of them. The gauge 
degrees of freedom reduce further (a 2 + b 2 + c 2 ) parameters. Then the dimension of 
the moduli space is 



dim c (A4' a 



b,c) 



a + c + lab — a 



(4.30) 



which implies 



:-i) 



dime CM' J 



:-i) 



b+c 



(4.31) 



What is interesting here is that, by changing the variables as 



q = xyz, 



Q 1 = -x, 



Q2 = y, 



(4.32) 



we can write (4.29) as 



Z, 



crystal 



n— n- 



q m Qi 



(4.33) 



To see this, let us consider a molten configuration of the crystal as in the left picture 



of figure 35, and split it into towers of atoms as in the right picture. Each such 



tower includes several squares composed of a blue triangle, a red triangle and a white 



rectangle. The relation (4.32) implies that each such square contributes q to the 



partition function. Furthermore, some of the towers have an additional blue triangle 
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with or without an additional red triangle on their top. According to (4.32), an 
additional blue (or red) triangle contribute —Q\ (or — Q 2 ) to the partition function. 
Thus, in general, we have the following three types of tower: 

1. A tower which has k(> 0) squares without any additional triangles on its top. 
This contributes q k to the partition function. 

2. A tower which has k(> 0) squares with an additional blue triangle on its top. 
This contributes —q k Qi to the partition function. 

3. A tower which has k(> 0) squares with an additional blue and a red triangle on 
its top. This contributes q k Q\Q 2 to the partition function. 

The whole contribution from a molten configuration p is the multiplication of contri- 
butions from all the towers in p. For example, p in figure 35 contributes q & Q\Q 2 . Note 



here that any molten configuration p is a collection of some towers, but the converse 
is not true. Namely, not all collections of towers give a molten configuration of the 
crystal. In particular, for fixed k, a single molten configuration p can include at most 
one second type of tower, while the first and the third types have no such restriction. 
This follows from our melting rule of the crystal. Thus, the first and the third types 
are "bosonic" while the second type is "fermionic." Summing up all the molten con- 
figurations is now equivalent to considering all the collections of the "bosonic" and 
"fermionic" towers. The bosonic towers contribute 



oo . oo _. 

TT^TT — - 



q m QiQ2 
to the partition function, while the fermionic towers contribute 



(4.34) 



\l(l -q m Qi). (4.35) 



m=0 



Then we find that the total partition function is given by (4.33) 



Now, let us show that the partition function (4.33) is consistent with the wall- 
crossing formula for the BPS index. We use the same method as in [SSI [59, 60J 25 We 
first blowup the suspended pinch point singularity so that the divisor T> wrapped by 
the D4-brane is topologically C 2 . The toric web-diagram after the blowup is shown in 



figure 36 The blowup parameters are regarded as the Kahler moduli of the geometry. 



3 See also [61] for a pictorial representation of the wall-crossing phenomena of D4-D2-D0 states. 
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Figure 36: A blowup of the suspended pinch point singularity. We have a D4-brane 
on T>. 



Varying them generically changes the BPS degeneracy of wrapped D-branes, which 
is called the "wall- crossing phenomenon." In particular, when we take a large radius 
limit of the blowup cycle, the BPS partition function becomes very simple [5E1 1591160] . 
In the large radius limit of the two-cycle flx, the geometry near the D4-brane is C 3 . 
The two-cycle {3% now transversally intersects with the divisor T> wrapped by the D4- 
brane. A perturbative string analysis then tells us that no D2-brane on 0i can form 
a BPS bound state with the D4-brane|_j The partition function of the BPS index is 
then given by the D4-D0 partition function on C 2 : 



Zn 



oo 1 

nr 2 

71=1 



T 



(4.36) 



where q denotes a Boltzmann weight for DO-charge. 

We now shrink fli and go back to the singular limit of the Calabi-Yau three- 
fold. Between the large radius and singular limits, we have an infinite number of 
wall-crossing phenomena. Any such wall-crossing is associated with an appearance or 
disappearance of BPS bound states in the spectrum. Suppose that a BPS bound state 



with charge T decays into two BPS states with charge Ti and V2_ The charge and 
energy conservations imply that such a decay is possible only if 



argZ(r 1 )=argZ(r 2 ) 



(4.37) 



where Z(T) is the central charge of BPS states with charge I\ Since the central 



charge depends on the Kahler moduli, the equation (4.37) can be solved by the moduli 



26 Note here any a'-correction to the BPS condition is suppressed in the large radius limit. 
27 Any decay channel involving a non-BPS state does not affect the BPS index. 
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parameters, which gives us a real codimension one subspace in the moduli space. Such 
a subspace is called the "wall of marginal stability." In our case, we are only interested 
in the charge 



r = V + Y^ M kPk ~ NdV, (4.38) 



fc=i 



where (3 k is the unit charge for D2-branes on the A;-th blowup two-cycle, — dV is the 



unit DO-charge, and M k ,N e Z. We take the basis of the two-cycles as in figure 36 



We use the same symbol T> in (4.38) to denote the charge for our D4-brane. Since our 
D4-brane is non-compact, its central charge is divergent. We therefore regularize it to 
write Z{T>) = ^A 2 e 2l ' p , where A — > oo should be taken in the final expression. The 
phase if expresses the "ratio" of the volume and B-field of the D4-brane and is fixed 
so that < ip < 7r/4 throughout our discussion. This regularization was first given 
in [62]. By taking a suitable parameterization of the Kahler moduli space, the central 



charge of the D2-D0 states can be written as Z((3 k ) = z k , Z(—dV) = l. 28 



The possible decay channels relevant for the wall-crossings are then [60J 

2 2 

fi = V + J2(M k - m k )f3 k -(N- n)dV, T 2 = ^ m k (3 k - ndV, (4.39) 



which implies 



fc=i fc=i 



2 



Z(T X ) ~ ^AV^, Z(T 2 ) = J2 m kZ k + n. (4.40) 

fe=i 

Here m k , n express the charges of the D2-D0 fragment. According to the Gopakumar- 
Vafa invariants [63], 164"] for the suspended pinch pointPjthe non- vanishing BPS indices 
of the D2-D0 states turn out to be 

fi(r 2 ) = 1 for (m 1 ,m 2 ,n) = (±l,0,n), (0,±l,n), (4.41) 

n(r 2 ) = -1 for (mi,m 2 ,n) = (±l,±l,n). (4.42) 

Thus, we have walls of marginal stability for (mi,m 2 ,n) = (±1,0, n), (0, ±l,n) and 



(±l,±l,n). The locations of the walls are specified by solving (4.37). In particular 



when we fix z 2 = 1/2 then the walls are drawn in the complex 24-plane as in figure 



37, where we denote each wall by W™ 1 '™ 2 . Here Imz 2 = means the two-cycle /3 2 



28 Here z k € C expresses the volume and the B-field of the A:-th two-cycle when Im z^ is large. 
29 For the explicit calculation of the invariants, see appendix D of [60 . 
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Figure 37: The walls of marginal stability in the complex 24-plane with z 2 = 1/2- All 
the walls are semi-infinite line from singular points (red dots). Along the red arrow, 
we move Im^i from 00 to 0, crossing the walls W^° and W^ 1 for all n > 0. 

shrinks to a point. On the other hand, Re z 2 7^ implies a non- vanishing B-field on /3 2 , 
which is necessary to avoid massless singularities associated with D2-branes wrapping 
/3 2 . Between the large radius and singular limits, we move \m.z\ from Imzi = +00 to 
Imzx = 0. To avoid massless singularities, we tune Re^i so that < Rezi,< 1/2. 



Then we cross the walls W^° and W^ 1 for all n > (figure 37). 

Since our Ti has no positive integer greater than one which divides out Ti, we can 
use the so-called "semi-primitive wall-crossing formula" [65] . The formula tells us that 
the partition function changes at each wall-crossing as 

(r 2 ,r>Q(r 2 ) 



Z-tZll 



-i)< r ' r VQrQ 



(4.43) 



where Q\ and Q2 are Boltzmann weights for D2-branes on j3\ and P2 respectively. The 
bracket (r,r 2 ) is the intersection product of the charges. To be more explicit, it is 
written as 

2 
(r, r 2 > = Y, m ^V, f3 k ) = mi, (4.44) 

fc=i 

where (T>, /3k) is equivalent to the intersection number of the divisor T> and the blowup 
cycle (3k- By taking into account all the wall-crossings between the large radius and 
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Figure 38: Left: The toric diagram of C 2 /Z3 x C. ^3 and p^ represent the lattice point 
(0, 3) and (1,0) for respectively. Middle: The brane tiling. Right: The quiver diagram 
for D2-D0 states. 



singular limits, we find that the partition function at the singular limit is written as 

2m-D2-i>o = Z, I I / •',.;% = I I — I I '—^- (4.45) 



z f\ l -« nQ i =TT — TT 1 ~ qmQl ~ 
H 1 - Q n QiQ2 t\ x - Q n ii 1 - <rQiQ, 



which agrees with (4.33). 



4.4 Orbifold C 2 /Z 



N 



X 



In this subsection, we study two dimensional statistical models for D4-D2-D0 states 
on C 2 /Zjv x C, where the D4-brane wrapping on a non-compact divisor. The toric 
diagram of C 2 /Zat x C is represented by the lattice points (0,0), (1,0) and (0, N) 
in Z 2 . We listed the toric diagram, brane tiling and quiver diagram associated with 



C 2 /Z 3 x C in the figure 38 



Orbifold partition 

The first example is the D4-brane wrapping on the C 2 /Z^v associated with the lattice 
point (1,0). The blow up geometry of C 2 /Zat is given by the Ajy-i-type ALE space. 
Since the intersection product of each blow up two-cycle and the four-cycle wrapped by 
the D4-brane vanishes, wall-crossing phenomena do not occur [60]. This implies that 
the BPS index of the D4-D2-D0 states is independent of the size of compact two-cycles. 
In the large radii limit of the two-cycles, the D4-D2-D0 states can be counted in terms 
of g-deformed Yang-Mills theory [661 EZ]- The partition function of the g-deformed 
Yang-Mills theory is related to the instanton partition function of the Vafa-Witten 
theory on the Ajv-i-type ALE space, which is given by the iu(iV) x-characters 
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On the other hand, in the small radius limit of the two-cycles, the D4-D2-D0 states are 
described by our melting crystal model. The absence of wall-crossings implies that the 
melting crystal model should also reproduce the character of su(iV). We will explicitly 
show this below. 

From the brane tiling, we find that the supersymmetric quantum mechanics on 
D2-D0 states is given by the well-known 7iAr_i-type quiver [1] with superpotential: 

JV-1 

W = J2 trMXi-iXi-i ~ XtXi). (4.46) 

i=0 

The boundary NS5-branes for the lattice point (1, 0) are the first and second ones 



in the middle picture of figure 38 They generally intersect with each other at iV 
different points in T 2 . The D4-node * is located at one of them. The choice is related 
to the holonomy of the gauge field at infinity on the D4-brane, which is classified by 
Tii(S 3 /Z N ) ~ Zjv [22]. Suppose that our D4-node * is located at the 2-th intersection 
point for some i = 0, • • • ,N — 1. There is a chiral multiplet $j localized at the point. 
Then the D4-brane induces an additional superpotential 

^flavor = tT(J$iI). (4.47) 

We here fix i = and consider the corresponding crystal melting model. The other 
choices of i are realized by shifting the labels of the quiver nodes k as k — > k + % (mod 
N) in the final expression. 

The dimer model and the perfect matching associated with the lattice point (1,0) 



are shown in the first and the third picture of the figure 39 The perfect matching 
contains all the adjoint chiral fields $&, and the constraint (3.6) implies that they have 
the vanishing vev's: 

$ fc = 0, (k = 0---N -1). (4.48) 

After eliminating $& from the massless spectrum, the quiver diagram can be depicted 



as in the left picture of figure [40} The non-trivial F-term conditions are now given by 

X N _ 1 X N _ 1 -X Q X o + IJ = 0, X k _ 1 X k _ 1 -X k Xk = (MO). (4.49) 

Note here that J = also follows on supersymmetric vacua. The two-dimensional 



crystal for the D4-D2-D0 states is now similar to that in the right picture of figure 28 
The only difference is that we here have N different types of box. The melting rule of 
the crystal now implies that 
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(0,0) 



(1,0) 



(0,3) 



Figure 39: The first: The bipartite graph of C 2 /Z 3 x C, which clearly admits an 
isoradial embedding. The second: The reference perfect matching. The third: The 
perfect matching associated with the lattice point p±. The forth: The perfect matching 
associated with the lattice point p 3 . 



• Each atom can be removed if and only if its left and lower edges are not attached 
to other atoms. 

Each molten crystal is regarded as a colored Young diagram, in which every box has 
a color. The color of a box atom is determined by the quiver node associated with 
it. Since we have N quiver nodes besides the D4-node, there are N different colors of 
box. We call a box "A;-colored box" if it is associated with the k-th quiver node. We 



show an example of molten crystal in figure [40j white, blue and red boxes represent 
color 0, 1, 2, respectively. 

We here emphasize that the molten configurations of the crystal correspond to the 
torus fixed point set (orbifold partition) of the ADHM moduli space on the ALE space. 
Counting the melting crystals is equivalent to evaluating the Euler characteristic of 
the ADHM moduli space, which is given by the level-one character of su(N). 

Here, the partition function of the melting crystal is given by 

oo JV-1 



Z, 



crystal 



e (-i)^w^n 



X 1 



(4.50) 



do,'"<2;v_i=0 



V(d) 



k=0 



The second summation is taken over the set V(d) of all molten configurations with dk 
fc-colored boxes for k — 0, • • • ,N — 1. We denote by x k the Boltzmann weight for the 
A;-colored box. 



To see that (4.50) is equivalent to the character of su(iV), we first evaluate the sign 
factor determined by the dimension of the moduli space. The chiral fields Xk,Xk,I 



^N-^ 



and J contain ^2 k=0 2dkdk-i + 2d degrees of freedom. The F-term conditions (4.49) 
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X n X 



-V, 






Figure 40: Left: The reduced quiver diagram of D4-D2-D0 states of C 2 /Z 3 x C Right: 
A molten configuration for D4-D2-D0 states. Each box has one of the N = 3 type 
color associated with the Z3 orbifolding. 



reduce J2 i=0 d\ parameters, and the gauge transformations further reduce the degrees 
of freedom by d\. Therefore, the dimension of the moduli space is 



N-l 



dimc(.M d -) = J2 ( 2 44-i - 2d 2 k ) + 2d . 



k=0 



This means that the sign factor is 

(_l)dimc(M,r) = 1 



(4.51) 



(4.52) 



We next evaluate the summation over V(d) in (4.50). We count the number of 



colored boxes in a molten configuration. It is shown in [26] that there exists a one- 
to-one correspondence between the set of the molten configurations U^V(d) and the 
set 



)| = nN,n = 0,1,2, ■■ \ , (4.53) 



\k = (k , ■■■ , fejv-i) eZ N \Ys k i = °} x { Y 

i 

with the identification 

1 JV-l N-l 

d * = 2 E h ) + E k * + n - ( 4 - 54 ) 

j=0 j=i 

Here \Y\ denotes the total number of boxes in a Young diagram Y . Then the partition 
function of the melting crystal factorizes into two parts [26J: 



N-l 



•'crystal 



iEft^Efi'^^ E (,„...„.,_,. 



(4.55) 



fc : £. fc l=0 *=0 «=0 Y:|Y|=n7V 

If we identify the Boltzmann weights of DO-charge and D2-charges as q = xo^i • • • xn-i 
and Qi = Xi,(i = 1,- ■ ■ , N — 1) respectively, the former factor becomes 



JV-l 






1 v iV-l ,2 , v iV-l , 



N-l 



\ "* T |" a n,i(ni-n i+1 ) Qfii 



(4.56) 
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Figure 41: Left: The reduced quiver diagram associated with the lattice point p 3 . 
Right: A molten configuration. 



\N-1 



Here we change the variables as n« = Y^j=i k{. On the other hand, the latter factor 



E E 



[Xi---X N 



oo .. 

11 (1 _ q m)N 
m=l v y y 



(4.57) 



n=0 Y:\Y\=nN 

corresponds to the generating function of the Euler characteristics of the Hilbert 
schemes of points on the A/v-i-type ALE space. Therefore we find that the parti- 
tion function of the melting crystal is written as 



Z, 



crystal 



JV 

q24 



77(g) 



A' 



J2 q^ nTCn Q n , 



(4.58) 



nez^- 1 



riV-1 



where C is the A^v-i-type Cartan matrix and Q n := Yl i=1 Q^ ■ This agrees with the 
level-one character of six(N) up to a Qfc-independent prefactorrj 

Another melting crystal 

We now turn to the second example, in which the D4-brane is wrapping on the non- 
compact divisor associated with the lattice point (0, A). The topology of the divisor 
is now C/Zat x C In order to be concrete, we will treat the example of N = 3. 
The lattice point corresponding to the divisor is denoted by p% in the left picture of 



figure [38} When we identify the boundary NS5-branes as the first and fifth ones in 
the second picture of figure 



38 



X 2 is now identified with Xp. Then the D4-node is 
now attached to the red and white quiver nodes, and induces the following additional 
superpotential: 



WWor — JX2I. 



(4.59) 



:in 



To be precise, we here obtain the level-one character for the trivial weight. The reason for this 



is that we set i = in (4.47). For general i in (4.47), we obtain the character for the i-th level-one 



weight of su(N), up to the n^-independent overall factor q a« e*-i *>i Vi with e Vj = Qj. 
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Figure 42: Left: A molten configuration of the crystal. Right: Division of the left 
picture into towers of atoms. A white, a red and a blue rectangle form a single square. 



The perfect matching associated with p% is shown in the rightmost picture of the figure 



39 Then the constraint (3.6) now requires 



X k = (A; = 0,1, 2) 



(4.60) 



on supersymmetric vacua. After removing X&, the quiver diagram can be depicted as 



in the left picture of 41 The non-trivial F-term conditions become 



$ 1 X -X % = 0, $ 2 X 1 -X 1 $ 1 = 0, %X 2 - X 2 $ 2 + / J = 0. 



(4.61) 



Again, the melting rule of the crystal is the same as for the young diagrams, but the 
coloring of boxes is different from the previous example; boxes with the same color sit 



in the horizontal direction (right picture in figure 41). 



The partition function of the melting crystal is given by 



Crystal = ^(-l)^"****- >Z* >Z* ^ 



d 2 

2 ■ 



(4.62) 



where p runs over all the molten configurations. The integers do, d\ and d 2 stand for 
the numbers of white, blue and red boxes in a molten configuration p, and Xo, X\ and 
x 2 represent the Boltzmann weights for white, blue and red boxes, respectively. The 
sign factor is determined as follows. The chiral fields contain the X/i=o(^*^»+i + df) + 



d + d 2 degrees of freedom. The F-term conditions (4.61) reduce ^2 i=0 did i+ i degrees 
of freedom, and further the gauge transformations reduce ^2 i=0 df degrees of freedom. 
Then the sign becomes 



f^\dim c (M d0}dl ,d 2 ) — (_\}<k+da 



(4.63) 



We can find the closed expression for the generating function by the similar manner 



to the second example in subsection 4.3 Let us consider a molten configuration as 
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in the left picture of figure |42l and divide it into towers of rectangles as in the right 



picture. We then change the Boltzmann weights as 

q = x xix 2 , Q\ = —xo, Q 2 = X\. (4.64) 

This means that a square composed of a white, a red and a blue rectangle contributes 
q to the partition function. An additional white or a blue rectangle contributes Q\ or 
Q 2 to the partition function. Thus we have the following three types of tower: 

1. A tower which has k(> 0) squares without any additional rectangle on its top. 
This contributes q k to the partition function. 

2. A tower which has k(> 0) squares with an additional white rectangle on its top. 
This contributes q k Qi to the partition function. 

3. A tower which has k(> 0) squares with an additional white and a blue rectangle 
on its top. This contributes q k Q\Q 2 to the partition function. 

The partition function can contain arbitrary number of the three-types of towers, and 

therefore can be written as 

00 -. 00 1 

Crystal = J] 7~^ II (1 _ q m Ql){l _ ^n QQ) ' ( 4 ' 65 ) 

n=l m=0 



To see that (4.65) reproduces the correct BPS partition function of the D4-D2-D0 



states, we consider the wall-crossing phenomena of BPS states. We consider the blow 



up geometry of C 2 /Z 3 x C as in figure 43 In the large radii limit of the two-cycles, 
the divisor wrapped by the D4-brane is isomorphic to C 2 . Then the D4-D2-D0 state 
counting reduces the D4-D0 state counting on C 2 , which again gives the partition 



function in (4.36). To obtain the partition function in the singular limit, we have to 
detect the positions of the walls of marginal stability and evaluate the jumps of the 
BPS indices. The possible decay channels are again written as the following form: 

r -> Ti + r 2 (4.66) 

with 

2 2 

T = V + J2M k p k -NdV, T 2 = J2^kPk + ndV. (4.67) 

fe=i fc=i 

The Gopakumar-Vafa invariants tell us that the non- vanishing BPS indices are now 
ft(r 2 ) = 1 for (mi, m 2 ,n) = (±1, 0, n), (0, ±1, n), (±1, ±1, n). (4.68) 
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Figure 43: The web diagram of blow up geometry of C 2 /Z,3 x C The grey region 
T> represents the toric divisor wrapped by the D4-brane. In the singular limit, the 
topology of T> becomes C/Z 3 x C The (5\ and /^ represent the blow up compact 
two-cycles. 



We have the walls of marginal stability for (1711,1712, n) 



The position of the walls are identified by solving (4.37). When we fix z 2 = \, 



±l,0,n),(0,±l,n),(±l,±l,n). 

the 



locations of the walls are the same as in figure 37 However, the intersection products 



between (3k and T> are different from those in subsection 4.3 



CP,/3i> 



<Z>,Aj) = 0, 



(4.69) 



Using the wall-crossing formula (4.43), we find that the BPS partition function at the 
orbifold limit is written as 



Z\ 



D4-D2-D0 



*nw 



1 



m=0 



q m Qi){l-q m QxQ2 



00 _. 00 

11 13^ 11 7i 



_ - q" 

n=l m=0 



- q m Qi)(l - q m QiQ, 



(4.70) 



which coincides with (4.65). 



We here comment on the relation to the instanton counting on C/Z^ x C As 
the orbifold partitions are in one-to-one correspondence with the torus fixed points of 
the moduli space of instantons on C 2 /Ztv, the molten crystals we are considering here 
correspond to the fixed points of the moduli space of instantons on C/Zjv x C [69J . 
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A Stability condition 

In [I], the slope function is defined by 

y\ c ^ tie dim Mi 
MM) := ^? q ° \. - , (A.l) 

and M is called ^-stable if every non-zero proper sub module N <Z M satisfy [i#(N) < 

jjl$(M). The ^-parameters are set in j3] so that #* = 1 and •&}. = for k e Q - We 

here describe that the ^-stability we use in this paper is equivalent to this ^-stability. 

First of all, the ^-stability is invariant under the following two types of change 

• d -» d + C for any (eM, 

• d — > £0 for any £ G M + . 

We particularly use the first one. For a given A-module M, let us define •&' :— $ + C 
with 

V, c p; i^dimM^ 
C " £< e& dimM, • ^ 

Then the ^'-stability of M is equivalent to the ^-stability of M. Note that we have 
H$i(M) = 0, and therefore //#(iV) < /x#(M) is equivalent to fi$>(N) < 0. 

We now identify our 6 as # = $'. Although our slope function 0(M) is different 
from /j,g(M), we can easily show that 9(N) < is equivalent to He(N) < 0. Hence, 
our ^-stability is equivalent to the original ^-stability. Furthermore, it follows from 
9 = d + C that 0* = 1, -& k = implies 0» > 0, fc < 0. 

B Isoradial dimer model 

The necessary and sufficient condition for a dimer model Q y on T 2 to admit an isoradial 
embedding was given in Theorem 5.1 of [70]. The condition is rephrased in terms of 
zig-zag paths as followsp] 

31 For the definition of zig-zag paths, see subsection 3.5 
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Figure 44: An example of dimer model which does not admit any isoradial embedding, 
which is associated with the local P 1 x P 1 . 



• Any zig-zag path of Q v is a closed curve without self-intersection, 

• In the universal cover Q v , any two zig-zag paths share at most one edge. 

All the dimer models shown in this paper admit an isoradial embedding. On the 



other hand, for example, the dimer model in figure [44] does not admit any isoradial 
embedding. In fact, there are pairs of zig-zag paths in the universal cover which share 
an infinite number of edges. 



C Perfect matching and toric divisor 

When we consider a single DO-brane probe on the Calabi-Yau singularity Ys, the 



condition (3.6) is understood as the condition that the DO-probe is moving in the 



divisor T>. We here describe this. 

Let us first note that the chiral multiplets X a e Q\ are expressed by commutative 
complex variables if we only have a single DO-probe. In this case, the F-term conditions 
without flavor branes are known to be solved by [HJ [321 E2] 



X, 



IJ Pm, 

m^X a 



(C.l) 



where the product runs over the perfect matchings involving X a , The variable p m is a 
complex variable and regarded as a new "field" associated with the perfect matching m. 
These new fields in fact trivialize all the F-term conditions dWo/dX a = 0. Therefore, 
we can regard p m as the fields of a gauged linear sigma model without superpotential. 
The moduli space of the sigma model corresponds to the background Calabi-Yau 
geometry Y% in which the DO-prove is moving. 
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Now, let us consider the condition (3.6). In the DO-probe setup, the condition is 
equivalent to 

Pmv = 0. (C.2) 



It was shown in [72 J that. 32 when T> is associated with a corner of the toric diagram, 
this condition describes the divisor T> in the moduli space of the gauged linear sigma 
model. This is physically interpreted to mean that the DO-probe should move in the 
divisor T> in order to keep the BPS condition with a D4-brane wrapping on T>. 
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